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ABSTRACT 


A  fundamental. problem  in  laser-beam  propagation  through  air  is  to  deter¬ 
mine  the  statistics  of  the  fluctuating  irradiance.  To  date  this  problem  is 
still  unsolved  in  the  strong-optics  regime  (where  scintillation  is  appreciable). 
This  work  presents  an  exhaustive  survey,  correction,  and  expansion  of  the  re¬ 
sults  that  can  be  obtained  through  selective  summation  of  averaged  products  of 
Bora-series  te/ms  substicuted  for  the  electric  field  in  the  average  field,  the 
mutual-coherence  factor,  and  the  N-th  moment  of  the  irradiance.  Detailed  error 
analyses  and  regimes  of  validity  are  provided. 
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GLOSSARY 


1.  COORDINATES,  VECTORS,  ETC. 

r  =  coordinate  or  spatial  vector. 
zm  =  propagAtion-direction  coordinate. 

P  “  two-d inens tonal  vector  in  a  plane  z  =  const. 
k  ■  Four ier-tr and form  conjugate  variable  of  r;  a  wavevector. 
t'a  =  Fourief-transfom  conjugate  variable  of  z. 

V 

K  -  Fourier-transform  conjugate  variable  of  p. 
r,  k  *=  unit  vectors. 

%  *■  transverse  derivative  x3/3x  +  y3/3y. 

=  transverse  Laplacian  32/3x2  +  32/3y2, 
d2r,  d2p  =  differential  area  in  z  =  const,  plane, 
d^r  =  differential  volume. 

67O  ■  two-dimensional  Dirac  delta  function. 

63O  =  three-dimensional  Dirac  delta  function. 

2.  ROMAN  ALPHABET 

A  -=  amplitude  of  electric  field  E. 
b  =  beam-wave  parameter,  see  (2.8). 

B  =  normalized  field  E/Eu. 

Cr2  =  refractive- index  structure  constant  (dim:  £"2/3). 

d  *  total  length  of  turbulent  medium. 

E  =  ilectric  field. 

E0  *  electric  field  in  free  space. 

r()  =  any  functional  dependence  needed  only  in  local  ;ontext. 

Fn  •=  exponential  factor  expressing  deviation  from  WKJB  form  of 
the  electric  field  (Section  IB)* 

G(rp,  r 2)  -  modified  Green's  function  or^spherical-wave  field  at  rp  due 
to  delta-function  source  at  ?£• 

G0(rp-r2)  =  free-apace  Green's  function. 

I  =  irradiance,  flux  or  any  quantity  « | E j  ^ . 

k  =  the  wavenamiier  ui/c  of  the  monochromatic  wave. 


L  =  pathlength. 

Lq  =  macroscale  of  turbulence. 

£0  =  microscale  of  turbulence. 

=  integral  scale  of  turbulence  [as  defined  in  (7.4)]. 
b(£)  =  number  of  eddies  of  size  i. 
n(£)  =  density  of  eddies  of  size  l. 
n(r,t)  *=  refractive  index. 

-a.  -a. 

Qm  =*  a  sum  of  Kj  from  j  =  m  to  j  =  n  in  the  n-th  Born  term  of  E. 
qm  =  a  similar  sum  of  kj . 

(t)  *  randomly  moving  center  of  eddy,  t  *  time. 

V  =  volume  of  turbulent  medium. 

3.  GREEK  ALPHAEET 

a  *  inverse  of  coherent-field  decay  length,  a  = 

A  *  Laplacian  operator,  cr  small  increment,  or  small  error. 

Abc  «  binary-cluster  error  ^  k^L^e^. 

A q  *  small-angle  scattering  error  'v  L  Km^^L0_^^c 

As  =  sagittal  error  ^  L2  Km1^3L0''8',,\_1e2. 

e(t, t)  =  relative  dielectric  permittivity  (dimensionless) . 

5e  =  deviation  of  e (? , t)  from  its  mean. 

e2  _  variance  of  e,  mean  square  of  6e. 

&  “  operator  notatio.i  for  6e. 

6e(K,z)  -  transverse  Fourier  transform  of  fc(p,z) 

n(x)  =  Fourier  transform  (three-dimensional)  of  6e(r). 

nj_(x)  =  Fourier  transform  (three-dimensional)  of  C^(r), 

•*  \ 

0  =  ai.gular  deflection  of  a  ray,  8  =  dp  /dz  (Section  8). 

<$0m  =  Km/k  =  scattering-angle  sine  at  rm  (Section  6). 

<m  =  5.92  £0-^  =  microscale  wavenumber. 

X  =  wavelength. 

5^(r)  =  shape  function  of  an  eddy  (Section  10). 

II  =  notation  for  a  product  of  factors. 

2  =  notation  for  a  sum  of  terms. 

=  a  variance  parameter,  an  expansion  parameter  (Section  4). 


aj2  = 
*00  - 

(j>  a 

X  - 

\J)  a 
* 


the  variance  of  x»  the  log-amplitude,  in  the  Rytov 
approximation  («  c2k?/6 

the  variance  of  I  normalized  to<I>  2. 

<5>  (K.k^)  *  normalized  spectrum  of  6e(t)* 

short-hand  notation  for  $(K,0). 

partial  spectrum  of  6e :  yields  $  (K,!^)  when  further 
transformed  with  respect  to  Az. 

the  phase  increment  of  E  above  kz. 

logarithm  of  the  amplitude  A. 

logarithm  of  the  field  E  (<{>  *  X  +  i  <f>) . 

the  phase  exponential  used  for  evaluating  points  of 
stationary  phase  {Section  6). 
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1 .  INTRODUCTION 


Laser-beam  propagation  through  turbulent  air  has  been  studied  seriously  for 
more  than  ten  years.  Nevertheless,  a  fundamental  problem  -  that  of  determin¬ 
ing  the  statistics  of  the  irradiance  of  even  the  most  elementary  light  wave 
(a  monochromatic  plane  wave  excited  at  z=0  and  travelling  to  z=L)  -  has  not 
been  solved.  The  difficulty  lies  specifically  in  the  region  of  parameters 
(Cn2,  L,  k)  where  amplitude  or  irradiance  scintillation  i3  appreciable.  Con¬ 
sequently,  standard  small-parameter  perturbation  techniques  fail  to  provide  a 
basis  for  calculation.  The  problem  is  so  very  fundamental  because  the  operation 
of  laser  beams  in  turbulent  air,  for  example,  is  limited  by  irradiance  fluctu¬ 
ations.  Laser  beams  can  be  considered  as  linear  superpositions  of  spherical 
waves,  hence  the  fundamental  problem  is  the  behavior  of  a  spherical  wave  in 
turbulent  air.  Nevertheless  spherical  waves  do  not  differ,  essentially,  from 
plane  waves  (except  at  the  origin  which  is  excluded  from  the  frbulent  region) , 
and  therefore  plane  waves,  which  are  somewhat  simpler  to  handle,  are  considered 
first.  To  our  knowledge,  three  basically  different  -  although  not  unrelated  - 
methods  do  form  a  basis  for  analysis  of  this  very  difficult  classical  field- 
theory  problem. 

(i)  The  development  of  approximate  equations  for  the  electric  field  E, 
or  for  related  quantities  such  as  InE. 

(ii)  The  development  of  equations  for  moments  of  the  electric  field,  e.g., 
for  the  mutual  coherence  factor,  the  four-point  correlation,  etc. 

(iii)  The  analysis  and  selective  summation  of  terms  obtained  by  substituting 
the  Born  series  for  each  field  E  in  any  moment  such  as<IN>  ,  the 
average  of  the  N-th  power  of  the  irradiance  I. 

Each  of  these  approaches  has  varyirg  developments  and  adherents  in  the  U.S.A. 
and  the  U.S.S.R.  Each  has  advantages  and  disadvantages.  Rather  than  making  a 
detailed  comparison  and  enumeration,  we  state  that  method  (iii)  is  the  one  we 
will  follow  throughout  this  work.  Its  one  great  advantage  over  the  other 
methods  is  that  it  allows  a  cumulative  error  analysis  of  the  simplifications 
and  approximations  made  on  the  way  to  results.  Thus,  it  allows  one  to  set 
bounds  or  regions  of  validity  for  the  solution  of  each  problem. 

This  work  will  contain  results  for 

(a)  the  coherent  (or  average)  field  <  E>, 

(b)  the  mutual-coherence  factor  <  E0?i)E*  (?2)>  , 

(c)  the  N-th  irradiance  moment  <  lN>, 

for  plane  waves  and  homogeneous  stationary  turbulence.  Some  of  these  results 
have  been  derived  previously  by  us  or  by  other  authors  using  different  methods. 
All  of  the  results  can  be  derived  for  spherical  waves  (with  appropriate  modifi¬ 
cations)  but  we  consider  that  a  non-essential  problem  here,  and  therefore  we 
restrict  all  derivations  to  plane  waves.  Lixewise,  many  applied-propagation 
results  are  of  no  interest  to  us  here.  Many  of  the  latter  parts  of  the 
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derivations,  e.g.  steps  leading  from  (7.2)  to  (7.4)  have  been  given  so  often 
and  in  so  nwiy  places  that  the  intermediate  steps  are  omitted.  The  motivation 
for  this  work  is  to  update  previous  work[7),  aj.so  reported  in  the  literature*, 
tc  include  the  optics  regime,  to  correct  errors  in  interpretation,  to  give 
extended  error  analyses,  and-  above  all-  to  give  an  exhaustive  and  comprehensive 
survey  of  the  selective- summation  method  in  computing  results  (a) ,  (b) ,  anA  (c) . 
Although  the  strong  optics  regime  where  k 7/b  1,11/6  Cn2  £  1  is  not  yet  understood 
as  far  as  irradiance  statistics  go,  the  selective-summation  method  appears 
promising  for  further  analysis. 

Ihe  organization  of  this  report  is  as  follows :  Sections  2  and  2  discuss 
various  forms  of  the  wave  equation  in  integral  form,  and  of  the  Born  series. 
Section  4  is  an  aside  in  which  the  cumulative  effects  of  small  errors  A  in  the 
Born  terms  are  discussed.  The  elimination  of  backscatter  (i.e.  of  radiation 
scattered  at  more  than  tt/2  radians)  is  the  subject  of  Section  5,  and  Section  6 
discusses  the  elimination  of  all  but  sm^ll-angle  scatterings  so  that  the  so- 
called  parabolic  or  bean-wave  Equations  (6-10)  and  (6-12)  ensue.  Geometrical- 
optics  approximations  yielding  the  WKR  expression  and  the  ray  equations  are 
detived  in  Sections  7  and  8.  A  most  important  result  of  Section  8  is  the 
discussion  of  the  variance  of  the  scattering  angle  6  in  multiple  integrals 

d^Km - d^Kn.  With  thi-3,  we  treat  the  wave  equation  statistically  in  Section 

9  to  prove  in  Sections  11  and  12  that  the  n-point  correlation < fie (1)  -  fie(n)> 

can  be  replaced  by  <5e(l)  5e(2)> - < fie (n-1) fie (n)>  in  all  moments  of  field 

and  irradiance.  We  develop  a  model  for  fie  in  Section  10  to  prove  this  most 
important  point.  The  mode]  can  be  tailored  to  fit  Kolmogorov  or  other  turbulence 
statistics.  Diagram- technique  bookkeeping  of  terms  of  the  Born-series  products 
are  introduced  in  Section  13,  the  average  field  <  E>  is  deduced  in  Section  14, 
and  the  mutual-coherence  factor  <  E(l)  E*(2)>  in  Section  15.  Sections  16-18 
are  devoted  to  higher-order  moments  <fN>  with  N  >1.  The  bookkeeping  of  terms 
is  explained  in  Section  16.  The  modif ied-Rytov  ai  "•roximation  and  its  range  of 
validity  are  deduced  for  the  irradiance  in  Section  17.  It  contains  energy 
•'onservation  automatically.  The  radiowave  case  is  discussed  in  Section  18,  and 
it  is  shown  that  a  Rice  distribution  for  I  holds  when  L  >>kLo^.  Finally,  the 
results  and  their  regimes  of  validity  are  summarized  in  Section  19  by  introduc¬ 
ing  a  two-dimensional  graph  of  the  permittivity  variance  e ^  vs  pathlength  L. 


*D.  A.  de  Wolf,  Radio  Sci .  3,  1379  (1967);'  3,  308  (1968);  J.  Opt.  See.  Am. 
58,  401  (1968);'  59,  1455  (1969). 
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2.  WAVE  EQUATION  IN  INTEGRAL  FORM 


Although  it  is  by  no  means  trivial  to  arrive  at  the  following  scalar  wave 
equation  for  high-frequenay  microwave,  and  optical  propagation  in  the  atmo¬ 
sphere,  we  assume  from  the  very  outset  that 

A  E  +  k2eCr,t)E  *  0  (2.1) 

where  E  is  the  electric  field  of  a  monochromatic  wave  propagating  through  air 
described  by  relative  dielectric  permittivity  e('f,t).  We  will  also  not  discuss 
further  the  conditions  under  which  dispersive  effects  are  properly  ruled  out, 
and  simply  set  e  =  n^  where  n  is  the  refractive  index.  Finally,  we  separate  off 
the  small  random  behavior  from  the  large  non-fluctuating  part  of  e,  i.e.  we  set 


e(r,t)  *  e (r'  [1  +  6e(r,t)]. 


(2.2) 


at  least  for  time  scales  in  which  e(r)  is  safely  assumed  to 
absorb  the  constant  (but  not  necessarily  uniform)  e(r)  into 
k.2  that  may  have  a  weak  and  slow  dependence  upon  r  compared 
Thus, 


be 

k2 

to 


constant,  and  we 
to  yield  a  new 
that  of  6e (r,t) . 


A  E  k2(l  6s)  E  »  0 


(2.3) 


is  a  convenient  shorthand  notation  for  the  scalar  wave  equation.  In  order  to 
deal  with  the  essentials  of  propagation  through  weak  random  media,  we  set  k  = 
constant,  and  restrict  ourselves  to  6e  <  1  everywhere  at  all  times.  Deviations 

from  k  =  constant  can  be  handled  but  they  add  greatly  to  the  complexity  of  the 
overall  iroblem. 

The  formula!  n  is  not  complete  without  specifying  initial  conditions.  We 
cnoose  the  form  one  out  of  many  possibilities)  that  field  E  reduces  continu¬ 
ously  to  E0  as  Se  ->  0  at  all  locations  and  that  E0  obeys 

2 

AE  +  k  E  =  source  term,  (2.4) 

o  o 

under  the  assumption  coat  the  source  term  is  very  restricted  in  space.  For 
example,  we  could  assume  a  plane-wave  source  at  r  =  (0,0,  -  00 ) ,  in  which  casa 
E0  “  exp(ikz).  Or,  one  can  assume  a  spherical--wave  source  at  the  origin,  in 
which  case  Ec  a  r“l  exp(ikr).  Although  we  have  set  the  righthand  side  of  (2.3/ 
equal  to  zero,  it  really  should  agree  with  that  of  (2.4).  However,  the  further 
formulation  will  be  for  E  -  E  so  that  the  source  term  drops  out,  and  the  equa¬ 
tion  for  E  -  E0  will  be  solved  in  terms  of  fie  and  E0  which  is  assumed  to  be 
fully  known. 
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Any  partial  differential  equation  such  as  (2.3)  plus  its  initial  and/or  bound¬ 
ary  condition  can  be  reformulated  as  an  integral  equation[l].  Standard  tech¬ 
niques  yield  as  equivalent  of  (2.3)  and  (2.4): 


2  /• 

Etr.t)  =  EQ(r)  J  d3^  GqG~v  )  Scfrj.t)  E(r,,t) 
G0(r"ri )  =  [exp (ik  !r-r1|)]/|r-r1| 


(2.5) 


This  equation  can  be  cast  into  a  form  of  ouch  greater  ue-  for  studying  whve 
propagation  in  a  random  medium,  namely  by  setting  E(r,t)  -  E0(r)B(r,  t) ,'  and 
attempting  to  find  an  equation  for  B(r,t),  the  normalized  field  that  describes 
the  effects  over  and  above  free-space  propagation.  This,  is  easily  done  by 
dividing  (2.5)  by  E0(r)  and  rearranging  the  factors  on  the  righthand  side.  The 
result  is 


B(r ,  t)  =  1  +  ^  5e(r1,t)  B(r1,t) 

G(r,r.)  =  G  (r-r.)  E  (r, ) /E  (r) 
i  o  1  o  1  o 


(2.6) 


This  equation  is  also  advantageous  because  B  is  dimensionless.  Furthermore, 
the  propagator  or  Green's  function  G(r,r^.i  takes  on  a  simple  form  for  plane 
and  spherical  waves.  Let  Gp  be  «-he  propagator  for  plane  waves,  and  Gs  likewise 
for  spherical  waves.  By  inserting  E0  ■  exp(ikz)  and  E0  =  r~l  exp(ikr),  respec¬ 
tively,  one  obtains 


G  (r,r  )  = 
P  1 


exp  |  ik  [|  r-rj^  |  -  (z-z^]) 


r  exp  |  ik[|  r-rx|  -  (r-r^]} 

r,r  )  =  - 1 - ~ - ^ 

1  t  Tb  I 


(2.7) 


rl  'r"rl 


Equations  (2.6)  and  (2.7)  constitute  the  integral  wave  equation  for  plane  and 
spherical  waves.  Likewise,  a  laser-beam  equation  could  be  produced  by  putting 
together  a  laser-beam  propagator  according  to  the  recipe  of  (2.6)  with[2]  the 
free-space  field  given  in  its  fundamental  mode  by 


E  =  -ikb  z-ikt 
o  \ 

-2  -2  -1 
b  f  r  -ik  R  , 
-  o 


^z-ikb^  exp  jikj^z  +  P^/2  (z-ikb^  | 


(2.8) 


where  p  is  the  vector  coordinate  in  a  plane  perpendicular  to  the  z  direction  of 
propagation,  i.e.,  r  =  (p,z),  with  rQ  the  beam  radius,  and  R  the  wavefront  radius 
of  curvature  at  z  =  0.  A  cursory  glance  at  (2.8)  suffices  to  convince  the  reader 
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that  the  laser-bean  propagator  is  more  complicated  than  those  of  (2.7).  How¬ 
ever,  it  is  well  known[2]  that  a  laser  beam  may  be  considered  as  a  linear 
superposition  of  spherical  waves,  and  therefore  the  fundamentals  of  the  prob¬ 
lem  of  laser -beam  propagation  can  be  understood  from  those  of  spherical  waves. 

Unfortunately,  the  spherical-wave  propagator  Gs  is  slightly  more  complex 
than  Gp  due  to  the  inclusion  of  some  additional  geometrical  factors.  So,  al¬ 
though  it  is  more  proper  from  the  point  of  view  of  laser  beams  to  work  with 
spherical  waves,  we  will  develop  the  formalism  for  plane  waves  because  that  is 
still  sufficiently  complex  to  justify  any  simplification  we  can  make.  Where  we 
specifically  require  spherical-wave  results,  we  will  make  the  necessary  modi¬ 
fications  in  reviewing  the  plane-wave  formalism  leading  up  to  the  result. 


3.  ITERATIVE  FORMS  OF  THE  WAVE  EQUATION 


The  integral  equation  can  be  expanded  into  a  Neumann-Liouville  series  by  itera¬ 
tion.  This  particular  N-L  series  is  also  known  as  the  Born  series  because  the 
generating  equation  is  in  essence  the  scalar  time-independent  Schrodinger  equa¬ 
tion  first  studied  in  this  form  by  Porn.  We  have  chosen  the  letter  B  for  the 
ratio  E/E0  for  that  reason.  Thus 

B  -£  Bn 

n=0  (3.1) 

n  ir  C  3  ■*  - 

B  =  n.  7-  Id  r  G(r  .,r  )  6e(r  ) 
n  m=l  4tr  J  m  m-1’  m  m 

with  the  specific  definitions,  B0  =  1,  r0  =  r  =  (p,  L),  and  with  the  appropri¬ 
ate  propagator  for  G(?m_i,rin) ,  Unless  otherwise  stated,  the  plane-wave  form  in 
the  first  of  (2.7)  is  implied. 

A  second  form  of  great  use  for  our  further  development  can  be  obtained  from 
(3.1)  by  utilizing  a  partial  Fourier  transform  of  6e(rm).  We  write  the  coor¬ 
dinate  rm  =  (Pn.» Zjn)  in  terms  of  the  Cartesian  coordinate  system  (x,y,z),  where 
z  is  the  direction  of  propagation,  and  indicate  the  (x,y)  coordinates  by  the 
two-dimensional  vector  p.  Strictly  speaking,  6e(r)  is  a  stochastic  variable, 
and  all  equations  are  more  generally  described  by  utilizing  Stieltjes  integrals, 
e.g.,  as  Brown[4]  does  in  his  work.  However,  in  order  to  avoid  a  more  formal 
treatment,  we  will  continue  to  use  ordinary  integrals  under  the  assumption  that 
6e( r,t)  varies  irregularly  but  smoothly  in  time.  Let  us  Fourier  transform  6e(r) 
with  respect  to  the  plane  coordinate  p,  and  define 

6e(K,z)  =  fd\  6e(p,z)  exp(iK.p)  (3.2a) 

The  inverse  transform  then  describes  the  synthesis  of  6e(_£)  in  any  plane 
z  =  const,  in  terms  of  two-dimensional  plane  waves  exp(-jK.p),  i.e., 

6s(r)  -  /“  K  ue(K.z)  exp(-iK.p),  (3  -  2b) 

4tt  *' 

with  a  weighting  factor  6e(K,z).  By  inserting  (3.2b)  into  (3.1),  we  obtain  the 
form, 


n 

b  =  n 


m=l  16tt' 


-  Id^r  fd^K  G(r  n,r  )  Se (K  ,z  )  exp(-iK  *p  )  (3.3) 

J  m  J  m  m-1’  m  m  m  r  m  m 


The  advantage  of  this  form  is  that  the  "scattering  problem"  can  be  solved  exactly. 
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By  this,  we  mean  that  we  consider  each  factor  of  (3.3)  as^a  three-dimensional 
integral  over  dz^  dZi^of  an  integral  over  d2pm  of  expO-iK^.p^)  GCr^,^) 
weighted  by  a  factor  The  d2pm  integral  represents  the  contribution 

at  r,^  of  the  electric  field  scattered  in  the  plane  z  -  zm  by  a  phase  screen 
expC-iKja.pju) .  The  6c(Km,zm)  integral  the..  adds  up  a  linear  superposition 

of  such  phase  screens  of  different  wavelengths  at  Zm,  and  the  dz,j,  integral  sums 
up  all  of  the  screens  as  Zj^  varies.  The  actual  scattering  part  of  the  five¬ 
dimensional  integral  in  (3.3)  is  thus  the  d^pm  integration.  The  rest  of  the 
problem  is  linear  superposition. 


Upon  inspection  of  G  from  (2.7)  -  we  use  Gp  from  here  on  -  we  note  that  it  is 
a  function  only  of  the  relative  coordinate  rm_i  -  rra.  Because  the  Jacobian  is 
unity,  we  may  replace  the  d2pm  integral  by  one  over  the  relative  transverse 
coordinates.  We  define 


m  m-l  m  m  '  i 

j=m 


Az  =  z  ,  -  z  , 
m  m-l  m 


and  utilize  the  transformation 


(3.4a) 


£  -  2;  [v?  -  v*J  • 


(3.4b) 


to  rtwrite  the  d  portion  of  each  factor  of  Bn  in  (3.3)  as 


2  r  0  exp  ik  |"(Ap  ^  +  Az  )  ^  -  Az  1  j 
I  ,2,  r  \  L  m _  nr _ m-jj 


k  i,2. 

71  Jd  A° 


2  .  2.1/2 
(Ap  +  Az  ) 
m  ra 


exp(iQ.n.Apm) 


(3.5) 


The  integral  can  be  performed,  see  Gradshteyn  and  Ryzhik[5]  formula  6.677, 
to  obtain 


(3.6) 


Consequently,  we  have  reversed  the  order  of  integration  in  (3.3)  and  we  are 
left  with  the  linear-superposition  integrals  (also  assuming  p  =  0) : 
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(3.7) 


oo 

Bn  *  fd\'  ■  'fd\  Jl  r!/dzm  dJ<Vz,,>  «  -  >!>2> 

8*  0 

x  exp  jik  [(l  -  ^/k2)172  |AzJ  -  Azj} 


-1/2 


This  is  the  basic  equivalent  to  Equation  (3.1).  We  shall  make  approximations 
from  here  on  by 'ignoring  parts  of  exponentials,  or  by  dropping  terms,  or  by 
modifying  factors  of  the  integrands  otherwise.  The  subject  of  approximations 
to  terms  of  a  series  is  one  that  is  very  essential  to  selective-summation 
methods,  and  the  entire  work  from  (3.7)  on  will  be  open  to  question  if  the 
approximations  are  not  justified.  That  is  to  say,  not  only  must  we  be  sure 
that  neglected  terms  are  small  compared  to  ones  retained,  but  also  that  the 
cumulative  effect  of  neglected  terms  is  small.  Therefore,  we  will  discuss 
error  analysis  again  as  the  formalism  is  progressed.  First,  we  consider  the 
general  cumulative  effect  of  dropping  or  si  (i>lifying  terms  in  the  Bom  series. 
An  outline  will  be  given  in  the  next  section. 
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4.  CRITERIA  FOR  APPROXIMATIONS  IN  THE  BORN  SERIES 


Consider  then  the  following  formal  representation  of  the  Born  series. 


B  =  £  (GS)n,  (4.1) 

n=0 

-v 

in  operator  language,  where  S  stands  for  6e(r)  in  the  case  of  (3,1),  or 
6e(K,z)  in  the  case  of  (3.7),  and  G  stands  for  the  propagator  and  the  integra¬ 
tion.  The  approximations,  we  wish  to  study  are  of  two  types: 

(I'i  Simplification  of  the  propagator'.  We  simplify  propagator  G  to  G. 
Then,  the  approximate  B  is  given  by  replacing 

(G  S)n  by  (G£)n  (4.2) 


(ii)  Simplification  of  the  statistics'.  Since  we  wish  to  compute  the 

average  field  <B>,  and  moments  of  the  irradiance  I  =  I0  BB*,  namely 

<  IN>  for  N  ■=  1,2  - ,  we  will  make  further  approximations  by 

replacing 

<  (G  6 ) 2n>  by  <  (G  6 ) 2>  n  (4.3) 


After  these  two  simplifications,  we  assume  that  a  closed-form  result  is  found 
by  rearranging  the  infinite  number  of  terms  in  <_B>  and  <  I^>  into  a  series 
of  powers  of  a  parameter  proportional  to  <(G£)^>,  and  performing  the 
summation. 


<B>  =  £  b  a 

n-0  °n 

<  1  >  ■  £  V,  ° 

n=0 


(4.4) 


In  other  words,  it  is  assumed  that  we  can  obtain  the  summable  results  (4.4)  by 
simplifying  G  to  G,  replacing  the  2n-point  correlation  <(Gg)in>  by  an  ordered 
product  (we  are  dealing  with  integral  operators!)  of  two-point  correlations, 
and  finally  by  rearranging  the  summation.  Consider  in  a  little  more  detail  an 
error  analysis  of  the  above  outlined  procedure  (the  decails  will  of  course  be 
described  in  subsequent  sections) . 

Ad  (i)  :  We  assume  that  the  propagator  G(0m,  A^)  of  (3.7)  is  simplified 
by  ignoring  a  small  part  of  the  exponent  and  of  the  denominator. 
Formally , 


where  A^  <<  1  and  A2  <<  1  are  small  variable  functions  of  Qm  and  Az^.  The 
error  of  approxin..  ting  G  by  G  will  be  estimated  by  choosing  an  upper  bound 
A*  to  all  Ai  +  A2,  and  by  defining  G(A')  =  G(l+A').  Because  A'  is  an  upper 
bound  and  a  constant,  the  factor  (1+A')  may  be  considered  to  be  a  numerical 
factor,  not  an  aerator.  Because  A'  <<  1,  the  error  in  <B>can  be  estimated 
by  comparing  <  -KA')>  to  <  B>  ,  and  likewise  for  <  lN>. 

Ad  (ii) :  It  is  more  difficult  to  explain  briefly  the  effect  of  statistical 
simplification.  The  approximation  (4.3)  amounts  to  the  neglect  of 
all  t  Tiulants  of  <(G£)2n>  higher  than  second  order,  l.e.,  for  a 
four- c-oint  correlation,  wa  neglect  the  effect  of  <(Gg)^>  - 

<  (G  *.  )2>2  (note  that  there  are  three  terms  contributing  to 

<  (G  ;  }->  2,  namely  so  that  S^S2>  £].£  3»  an^  £4  are  paired  to- 
geth  'i.  In  ignoring  higher-order  cumulants,  we  are  dropping  terms. 
The  i-tor  in  so  doing  is  assumed  to  be  (proof  will  be  given  in 
Sect!  ■  7  11  and  12), 


<  (G  fi )  2>  n  [l  +  0(A3)]n. 


(4.G) 


Consequently  wt  sec  that  the  error  can  be  estimated  by  finding  an  upper 
bound  o'  to  all  and  investigating  the  difference  between  results  obtained 
by  jipproxiraatiri  >  (Gg)2n>  first  by  <(Gg)2>n  and  then  by 
<(G£)2  (It  A"  >\ 

Thetwo  approxi  .1  ions  (i)  and  (ii)  are  then  summarized  by  replacing  G  first 
by  G,  and  then  n;-  5(A)  =  G(1  +  A),  where  A  =  A’  +  A".  Because  <  (G&)2> 
corresponds  to  t  *  in  the  rearrangement,  <(G8)2  (1  +  A)2>  corresponds  to 
a2  (1  +  A)2.  Tb  \ afore  the  error  in  <B>  is  estimated  by  comparing  the  fol¬ 
lowing  two  serie  r 


<  B  >  ~  5^  b  a2n 

n=0  °n 

<sw>  -  E  <>„„  [o(i  +  4>]2n 

n=0 


2-2 

Thus  if  <  B>  is  tummable,  it  is  a  function  of  a  ,  namely  <  B>  =  f  (0  ).  The 
second  of  Equations  (4.7)  then  states  that  <  B(A)>  =  f  [o2  (1  +  A)2].  The  error 
in  <  B >  is  then  estimated  by 


r~ 


02fl 


-»  A) 


a2[l 


+  G'A)] , 


(4.8) 


.  i.e.  by  regarding  the  effect  of  terms  of  0(A)  added  to  the  parameter  o2 
wherever  it  occurs  in  the  final  approximately-summed  expression.  The  same 
analysis  holds  for < IN>  .  Thus,  the  condition  A  <<  1  is  sufficient,  even  in 
a  cumulative  sense,  to  guarantee  that  we  car  utilize  x  =  in  f(x)  rather  than 
o2(1-A}2  <  x  <  o2(l+A)2.  Later,  we  will  identify  A  with  nL0/L  and  with 
n<(66)2>  ,  i.e.  with  errors  proportional  to  the  running  index  n  in  (4.7).  The 
procedure  will  be  modified  to  deal  with  this  c£se  ir  Section  6. 


5.  ELIMINATION  OF  BACKS CATTER 


From  this  point  on,  wa  shall  be  interested  in  random  madia  in  which  the  wave- 
numbers  K,  which  define  the  Fourier  decomposition  6e(K,.,)  of  the  permittivity 
variation  in  the  plane  perpendicular  to  that  of  propagation,  are  essentially 
restricted  to  2ir  L0"l  <  K  <  2ir£0“l  where  £0  and  LQ  are  the  micro-and  macro¬ 
scales  of  turbulence,  and  2tt£0~1  <<  k.  In  optics,  the  latter  inequal¬ 
ity  is  well  satisfied,  but  in  radiowave  propagation  it  is  m.  re  often  the  case 
that  k  £  2it£0-1.  In  both  cases,  however,  2ttL0-^  <<  k.  Consider  for  the  present 
that  2it£0“1  <<  k.  Then,  a  first  glance  at  (3.7)  indicates  that  the  factor 
(1  -  Qjn2/k2'i  is  quite  close  to  unity,  except  for  the  pathological  ceses  in  which  m 
is  very  large,  and  the  directions  of  all  Kj  for  j=m  to  j=n  (small  m)  are  such 
chat  the  small  vectors  add  up  constructively  to  exceed  k  in  length.  When 
(1  -  w  1,  the  last  factor  of  (3.7)  will  oscillate  slowly  for  A%  >  0, 

namely  over  a  length  scale  ^  k/Qm^,  and  rapidly  for  Azm  *  0,  namely  ever  a 
length  scale  k“l.  In  the  latter  case  of  Azra  <  0,  the  propagator  oscillates 
sinusoidally  many  times  over  a  length  scale  in  which  dt^m.Zjn)  changes  with  % 
(namely  a  distance  at  least  of  order  £0  and  more  probably  of  order  £:  an  appro¬ 
priate  mean  scale  between  £0  and  L0  such  as  the  integral  scale[6]).  The  oscil¬ 
lations  tend  to  cancel  the  contribution,  and  it  therefore  appears  reasonable  to 
restrict  the  dzm  integration  to  the  forward-scatter  regime  0  <  Zrn  <  Zm-1>  The 
lower  bound  of  the  dzjn  integration  is  zero  because  we  set  <5e  *  0  for  z  <  0. 
otherwise  we  lave  a  non-physical  situation  in  which  electromagnetic  energy  has 
propagated  an  Infinite  distance  before  reaching  z  (a  violation  of  all  the  approx¬ 
imations  we  '.nail  make).  Similarly,  we  should  set  6e  =  0  for  z  >  d,  the  medium 
width. 


Analytically,  this  is  expressed  by  taking  under  consideration  the  dzm  integra¬ 
tion  of  all  factors  in  Bn  dependent  on  z^.  This  part  of  the  integrals  can  be 
written. 


m-1 


-ik,.  tz 


f  dzm  6e(Km,zJe  f  m 


m  mm 


/n  -ik,  a; 

dz  6c(K  ,z  )e 
m  mm 


=  k 


^  k 


1  -  (1  -  Q  2/k2)1/2 
;  m 

1  +  (1  -  Q  2/k2)1/2 
m 


m 


m-1 


(5.1) 


There  is  a  factor  containing  %  in  the  (m+l)-st  factor  of  Bn,  but  we  need  not 
consider  it  here.  Each  term  of  (5.1)  is  Gaussian  to  goou  approximation  with 
zero  mean  and  a  variance  that  we  compute  by  ignoring  small  boundary  effects 
and  utilizing  the  properties. 


<fit(Ki,z1)  6e(K2,z2)> 


Mk'Az 

J>2  (i’ ,  Az)e 


c  '  -1>2(K},az)  6,/k,  +  K,) 

=  'KK.O  =  i  [(K2  +  k,2l1/2]. 


(5.2) 


— .00 
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The  first  term  of  (5  1)  has  a  variance  proportional  to  *(Kfn,kf),  and  the 
second  term  has  one  proportional  to  (d  -  z^j)  $0^,^) .  As  each  term  is 
Gaussian  with  zero  mean,  the  contribution  of  each  is  not  likely  to  exceed 
several  times  its  variance. 

Wavenumber  kj,  is  always  greater  than  ki  and  when  Qm  <<  k,  it  is  close  to  2k. 
Therefore  the  argument  of  $(1^,^)  -  which  is  (Km2  +  ki,2)l/2  —  always  exceeds 
k,  and  when  Qq  <<  k,  it  is  close  to  2k.  In  optics  it  is  thus  quite  permissible 
to  drop  the  backs  matter  term  because  there  are  essentially  no  components  to  •> 
at  wavenumbers  so  high.  The  very  small  corrections  to  zero  are  of  the  cruer 
&2  «  1  described  in  (4.5),  hence  they  do  not  amount  to  anything  appreciable. 

In  radiowave  propagation,  the  argument  for  ignoring  backscatter  is  slightly 
more  complicated.  It  goes  as  follows.  Consider  (4.4).  The  results  obtained 
for  radiowaves  all  indicate  (see  later)  that 

CO 

a2  *  J  dKK$(K)5  (5.3) 

o 


I*-  can  be  shown,  e.g.,  for  the  Kolmogorov  spectrum,  that  there  is  an  Z  with 
jtQ  <<  £  <  L0  such  that  the  above  integral  does  not  vary  appreciably  if  one 
sets  <KK)  -  0  for  K  >  2iri,“l.  However,  in  radiowave  propagation  ie  may  assume 
that  k  »  2irS.~l,  even  though  k  £  2iri0“^.  Since  none  of  the  remits  change 
appreciably  by  this  cut-off  (the  arguments  of  Section  4  include  cut-off  errors 
of  this  type)  we  make  it,  and  thus  have  set  $(Km,kjj)  *  0.  Therefore  backscat¬ 
ter  is  unimportant  in  this  case  too. 


as  a  good  approximation  for  .\  <<  Z,  where  Z  ■  lQ  in  optics,  and  Z  *  Z^t  the 
integral  scale  in  radiowave  propagation. 
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6.  SMALL-ANGLE  SCATTERING  APPROXIMATIONS 


So  far,  backscatter  has  been  eliminated,  and  the  Born  series  has  been  reduced 
somewhat  in  complexity  to  (5.4).  Further  reduction  appears  possible  in  view 
of  the  fact  that  when  index  n  is  not  too  large,  the  factor  (1  -  Q^/k2)  does 
not  differ  greatly  from  unity  because  K/k  is  small  for  most  (if  not  all) 
eligible  K-vectors.  However,  let  us  look  into  the  matter  more  analytically. 

We  reconsider  the  exact  derivation  of  (3.6)  and  (3.7)  from  (3.5)  by  an  approxi¬ 
mate  method  which  lays  bare  the  physical  background:  the  method  of  stationary 
phase  (one  of  the  well-known  methods  of  asymptotic  analysis) .  Reconsider 
(3.5)  in  the  following  form 


exp(iy 


+  Az 


m 


2jl/2 


m 


k  [ (Ap /  +  Az  ) 


1/2 


m 


m 


Az  1 
m 


(6.1) 


The  method  of  stationary  phase  states  that  the  integrand  contributes  appre¬ 
ciably  to  the  integral  only  in  the  direct  vicinity  of  the  stationary  points 
of  1%,  i.j. ,  the  locations  where  the  first  derivatives  of  Tm  with  respect  to 
A^  aiJ  Aym  are  zero.  Let  us  label  the  positive  square  root  of  (Ap^2  +  Az^) 

by  Ar  .  It  follows  from  the  above  that  the  stationar\  ooints  are  at 
m 


Ar 


m 


(6.2) 


It  is  then  easily  seen  by  further  application  of  the  conventional  stationary- 
phase  method  that  (3.6)  -  hence  (3.7)  -  follows  from  (6.1)  except  for  a 
slight  difference  in  the  denominator.  The  important  point  here  is  the  in¬ 
terpretation  of  (6.2).  Tor  given  Qm/k,  a  small  vector  let  us  say,  and  for 
given  pa_^  and  zm_^,  there  are  two  locations  (pm,zm)  for  the  m-th  scattering 
(i.e.,  for  the  m-th  stationary  points).  They  are  illustrated  in  Figure  1(a): 
one  is  for  a  "forward"  scatter  and  the  other  for  a  "back"  scatter  (^  •»  zm_^) . 
There  is  thus  a  straightforward  geometrical  interpretation  for  the  K-vectors. 
Let  us  name  6ra  3  -  Pm)/|r’m_i  -  fm|,  0m  is  a  vector  of  which  the  two 

cartesian  components  indicate  the  direction  sines  of  the  angles  0^  and  6my 
which  the  vector  rm,-i  -  makes  with  the  direction  of  propagation. 

The  discussion  in  Section  5  has  led  to  elimination  of  the  back-scatter  sta¬ 
tionary  po'nt,  the  one  for  whicii  the  direction  angle  is  larger  than  tt/2. _ ^ 

Only  forward-scatter  stationary  points  remain.  Consider  the  difference  60m 
of  •"he  direction  sines  atj^  aiu*  previous  point  r^.  ,  o0m  =  0m  -  0^^. 
From  this  definition  i.f  60m  and  from  (6.2)  it  follows  that 
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(b) 


Figure  1.  (a)  Forward-arid  back-gcatter  stationary  points  r  . 

(b)  Scattering  angle  60^  =  -  K^/k.  l' 


60  3  -  K  /k  '  r  K  -  -  k  60  (6.3) 

m  m  m  m 

•V 

The  interpretation  of  Km  is  analogous  to  that  of  momentum  transfer  in  quantum 
mechanics.  Consider  a  vector  Itm  of  length  k  in  direction  rm_i  -  rm,  and  a 
vector  km.  i  of  length  k  in  direction  rm  -  r^^.  By  (6.3)^the  vector  ^ 
has  a  component  perpendicular^to  the  z-direction,  namely  =  -  k66m.  In  the 
SchrBdinger  equation,  vector  represents  momentum  at  the  m-th  point  of 
scattering,  hence  is  the  momentum  transfer  perpendicular  to  propagation. 
When  6C«n  <<  tt/2 ,  it  is  approximFted  quite  well  by  the  actual  scattering  angle 
depicted  in  Figure  1(b)  (sin  60  555  50  for  small  angles),  with  error  of  the 
order  of  (<$Sm)  • 

For  optical  propagation,  the  inequalities  %  2vl0  ^  <<  k  are  well  satisfied. 
For  radio  propagation,  that  is  not  the  case,  but  the  outlined  err  analysis 
-uiL  Section  4  indicates  that  we  need  consider  only  those  for  which  Kj,  £  2irl 
wneic  "Z  <<  £  <  L0.  One  may  still  consider  <<  k  for  all  eligible  wave- 
vectors,  be-<iuse  2tt£-1  <<  k  in  this  case. 
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Now,  the  small-angle  scattering  npproximatioa  in  (5.4)  amounts  Vo  replacing 
(1  -  Qm^/k2)^/2  by  i  j[n  the  denominator  and  by  1  •  Qm2/2k2  in  the  numerator 
of  the  last  factor.  A  careful^discussion  of  "b-  ’.fications  of  this  re¬ 
placement  is  required  because  Qm  concains  n-m  va.  vectors  which  can  at!  d  up, 
conceivably,  to  a  length  k  if  n  becomes  large,  even  though  Kj  <<  k  for  all 
eligible  values  of  K j .  Such  a  discussion  does  not  appear  to  be  available 
elsewhere  at  the  time  of  writing,  except  -  in  rudimentary  form  -  in  an  earlier 
version  of  part  of  this  work[7]. 

It  is  helpful  to  return  to  the  definition  of  Qm  (for  fixed  n)  and  to  invoke 
(5.3),  so  that  is  written  as 

n 

Q^/k  =  Yj  SQj  (6.4a) 

j=m 


Although  it  is  difficult  to  define  exactly  what  is  meant  by  ..he  "behavior" 
of  Qm/k,  we  will  attempt  to  describe  it  by  noting  that  for  fixed  zm,  — , 

zn,  a  .function  f(Qm/k)  occurs  in  the  multiple  integral, 


u 


d260  6o(k60  ,  z  ) 

n  mm 


6e (k60  ,  a  )  f(Q/k) 
p  n  m 


(6.5) 


The  multiple  integral  is  random,  hence  we  are  interested  in  its  ^statistics. 


When 


Zj-Zj+i  exceeds  several  times  L0,  we  note  from  (5.2)  that  6e(W60j>  zj) 


and  <5e(k60j+i,  zj+])  are  uncorrelated .  It  ve  set  the  condition  nL0  <<  L,  as 
we  will  in  Sections  11  and  12,  it  then  follows  that  each  of  the  d2o0j  6e(k60j,  z^) 
integrals  is  uncorrelated  with  any  of  the  others  except  /or  a  fraction  *v 

(n-m)L0/L  of  the  n-m  dimensional  volume  defined  by  the  integrals  dzm  -  dzn 

over  the  fixed  locations  Zy  Under  tho  assumption  (proven  latei.  in  Sections 
11  and  12)  that  the  eifects  of  this  smalt  volume  are  negligible,  it  follows 

that  60m  -  60n  are  statistically  uncorrelated  in  the  sense  that  each  is 

weighted  by  an  independent  random  factor  6£(k60j,  Zj)  in  (6.5).  It  is  in  this 
sense  that  we  regard  Qm/k  as  the  sum  of  n-m  independent  random  variables  60 j, 
as  in  (6.4a).  Later,  we  shall  define  more  precisely  what  we  mean  by  this 
random  behavior,  but  now  we  note  that  Qm/k  can  therefore  be  considered  to 
approacn  a  Gaussian  random  variable  when  n-m  10  (for  smaller  values  of  n  we 
are  not  much  interested  in  the  cumulative  effects  an>way) .  Thus 


«C/k  > 


<  60. 

l 


66. > 

J 


=  (n-m)  <(60)  >. 


(6.4b) 


although  we  cannot  yet  specify  in  detail  what  the  averaging  brackets  specifi¬ 
cally  mean  in  this  formula.  Because  Qra  is  Gaussian,  when  n-m  not  small,  it 
follows  that  the  magnitude  of  Qra2/k2  is  restricted  to  seve.-.l  times  the  variance 
(6.4b).  Hence,  if  we  can  specify  more  preci.  1 y  what  the  averaging  procedure 
in  (6.4b)  is,  we  can  then  estimate  the  cumulative  magnitude  of  Qffi2/k2.  At 
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this  time  we  note  that  we  may  make  appt oximations  in  (5.4)  based  upon  utilizing 
«  1  under  two  conditions,  A  «  L 


(1)  A  »  nL0/L 
(ii)  A  =  n<(60)2> 


(6.6) 


2 

I.e.,  we  ignore  terms  of  0(nLo/L)  and  of  0(n  <(66)  >)  in  (5.4).  However,  we 
have  to  modify  the  discussion  of  Section  4  a  little  because  the  error  term  per 
factor,  A,  is  now  proportional  to  n,  whereas  the  analysis  assumed  A  independent 
of  n.  One  additional  assumption  will  amend  the  situation.  To  sea  this,  re¬ 
consider  (4.7)  and  define 


f(a2}  =  E  b  c2n 

n 

n=o 

N 

V‘2> 5  E 

n-o 


(6.7) 


2  2 

Obviously,  the  limit  of  fjg(a  )  is  f(a  )  as  we  allow  N  -*■  «.  However  the  new 
assumption  is  that  for  all  N  lar'*  jr  than  some  Nq  >>  1  we  have 


f(o2)  -  fN(o2)  'r  0(o2/N).  (6.8) 

2  2 

Obviously  this  holds  true  for  exponential  series  yielding  f(o  )  =  exp(+  a  ), 
and  it  will  be  useful  for  this  work  because  all  of  the  statistical  results 
are  simple  functions  of  such  exponentials, 

2 

Assumption  (6.8)  states  that  we  need  not  consider  any  Bn  with  n  »  o  for  the 
applications  in  mind,  i.e.,  fo*  calculating  <B>  and  <1^ > .  The  assumption 
Q^2  <<  k2  therefore  determined  by 

(i)  L  o2  «  L 
o 

(ii)  <(66)2 >cr2  <■<  1  a2  «  (kl0)2  optics  (6.9) 

2  2 

a  <<  (ki >  radiowaves, 

where  we  have  set  K  <  2u£.0  ^  for  optics,  and  K  <  2tt i?  ^  (i  ‘u  ,  the  integral 
scale)  for  radiowaves.  Thus  as  long  as  the  inequalities  (6.9)  are  satisfied, 
it  appears  to  us  valid  to  assume  /k2  <■  1  in  further  development  of  (5.4). 
Accordingly,  we  drop  Qm2/k2  in  the  denominator,  and  expend  (1  -  Qm2/k2)l/2  =• 

1  -  Qm2/2k2  in  the  numerator  to  obtain, 
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B 

n 


fd\  - 


zm— 1 


f*\  n  a 

*  m=l  8tt 


/  dz» 


6e(K,z) 
m  m 


“'K'tvi  -  *»>H 


(6.10) 


There  is  one  more  approximation  hidden  in  (6.10).  In  approximating  the 

numerator  of  (5.4),  we  have  to  make  sure  that  the  maximum  value  of  the  first 

ignored  term  in  the  exponent  does  not  add  up  to  an  appreciable  phase  effect, 

i.e.,  that  (Qm4  -  Ql+i^/Sk3  <<  1  t;r  each  If  we  set  Qm  =  +  Qnrt-1, 

then  we  obtain  the  condition  (using  z  <  L) 

m 


(K  +  4 


m 


K‘  K  *Q  .. 
m+1  to  m+1 


+  6  \2(Li +  4 


vw 


m+l 


)  «  1 


(6.11) 


The  inequality  K^L/k3  <<  1  can  be  written  as  L  <<  k  l  (optics)  and  L  «  k  i. 
(radiowaves) :  it  is  known  as  the  sagittal  approximation  in  optics  when  terms 
of  order  L/k3f^  are  dropped.  The  other  terms  represent  possible  cumulative 
effects  which  we  discuss  in  Section  14. 


Finally,  it  is  important  to  point  out  that  (6.10)  also  results  from  the 
parabolic  equation  favored  by  Russian  authors [8].  The  parabolic  equation 
follows  from  (2.3),  which  we  rew 'ite  as 


AB  +  2ik«VB  +  k  6e  -  0 


(6.12) 


by  substituting  E  =  B0  exp  ikz  (k  =  kz  =  kz/z).  This  equation  is  quite 
general,  and  equivalent  to  (3.7).  The  parabolic  aquation  is  obtained  from 
(6.12)  by  ignoring  o2g/3z2  £n  the  first  term.  It  is  equivalent  to  (6.10). 
However,  we  believe  that  the  discussions  of  the  parabolic  equation  in  the 
literature  (see  Reference  [8]  for  example)  do  not  make  the  nature  of  the 
approximations  sufficiently  clear,  hence  the  digressions  in  Sections  4-6  on 
this  topic.  We  will  return  again  in  Sections  8  and  14  to  the  difficult  topic 
concerning  the  nature  of  error  terms  implied  by  (6.9)  and  (6.11).  From  this 
point  on,  Equation  (6.10)  is  the  basic  starting  point. 
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7.  THE  WKB  OR  MOLIERE  APPROXIMATION 


The  Born  series  Bq  in  the  small-angle  scattering  approximation  (6.10) 
summable  when  the  exponential  function  is  replaced  by  unity.  In  that 
utilize  the  multiple- integral  properly 


xs 

case  we 


n 

n 

m=l 


m-1  r  L 

/  dzm  f(V  *  £f|  / 


dz  f(z) 


n 


(7.1) 


to  obtain  a  special  form  of  the  Wentzel-Kramers-Brillouin  (WKB)  approximate 
solution  of  a  wave  equation. 


B  =  exp 


=  exp 


Lt 

^2  j  dz  Jd2K  5e(K,z) 


o 

L 


rf 

L.  o 


dz  6e(p,z) 


p=o 


(7.2) 


also  known  as  the  Moliere  approximation [7 ,9] .  The  restriction  p=0  is  not 
essential.  The  errors  are  given  by  A  =  exp[-i  Qm2  (zm_i  -  zm)/2k]  -  1 
^  kL<(66)2>for  small  n.  Hence  the  error  analysis  requires 


a2  kL  <(69) 2  >  «  1 


(7.3) 


Before  further  discussing  the  error,  let  us  analyze  the  meaning  of  (7.2)  for 
a  random  medium.  It  states  that  the  electric  field  E  is  modified  from  E0  by 
random  phase  factor  exp(i<)>)  where  i<J>  is  the  exponent  of  (7.2).  By  virtue  of 
a  variant  of  the  central-limit  theorem  of  statistics  (which  is  valid  because 
L0  <<  L) ,  <j>  is  Gaussian  (with  zero  mean)  ,  hence  the  statistics  of  B  are  given 
by  the  variance  <<j>2  >  :  in  this  case  only<B>is  of  interest.  The  calculation 
of  <<f>2>  f  and  of  other  covariance  statistics  of  B  based  on  (7.2)  are  too  well 
known  to  repeat  here  (see  Tatarski's  work[10,ll],  or  any  of  the  more  compre¬ 
hensive  review  articles  on  wave  propagation  in  random  media) .  The  results 
are  (for  a  statistically  uniform  medium) : 


<B>=  exp[-  <4>  >/  2] 

co 

/ 


2  2 

,2  k  e  L 
<*  >_  Ttt 


dK  K$(K)  =•  2k?'i.L  e2 


<B(y  p.L)  B*(-  \  P,L)>  =  exp[-  \  D  (p)J 
D  (p)  =  <[<(>(|  p,l)  -  H-  \  p >l) ] 2 > 


i2  2T 

k  e  L 

4it 


/ 


dK  K<5> (K)  [1  -  I  (Kp)  ] 


(7.4) 
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with  the  usual  NBS-Handbook[13]  notation  for  Bessel  functions:  J0(Kp)  in  this 
case.  Extensions  for  media  in  which  mean  and  variance  vary  slowly  are  ob¬ 
tained  by  replacing  L  by  an  integral  of  dz  from  zQ  =  0  to  z  =  L  and  allowing 
e2  and  4>  to  depend  on  z. 

Error  Analysis:  The  inequality  (7.3)  is  more  stringent  than  the  second  of 
(6.9),  and  also  more  jo  than  (6.11)  -  as  will  be  shown  in  (14.8).  Let  us  there¬ 
fore  consider  it  the  defining  restriction  of  the  Moliere  approximation.  For 
£  1,  we  utilize  (7.3)  to  make  an  estimate  of  the  error  by  setting 
<(60)2 >£  (ki)-2  with  *  =  Thus, 

L  «  kL02  (7.5) 

2 

is  a  weak-effect  restriction.  In  order  to  obtain  the  restriction  for  a  >>  1, 
we  need  to  define  <(60)2><j2  in  (7.3).  We  do  this  in  the  next  section. 
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8.  BAY  EQUATIONS 


The  considerations  of  Sections  6  and  7  c?n  be  made  clearer  with  the  concept 
of  rays.  Usually,  the  ray  concept  is  treated  in  terms  of  a  differential- 
equation  formulation  of  the  wave  equation.  However,  we  can  also  do  it  in  the 
context  of^our  integral  formulation.  To  do  so,  we  return  to  (2.6)  but  rewrite 
B(*l)  =  A(r^)  exp[i$(ri)].  The  time  dependence  is  still  assumed  but  not  de¬ 
noted  explicitly.  Consequently, 


B  =  1  + 


a 

j  dz1  Jd2 p1  G(r  -  rx)  6e(r1>  A(r^) 


i<Kr1) 


(8.1) 


As  in  Section  6,  we  apply  the  method  of  stationary  phase,  this  time  to  the 
d^p^  integration,  and  the  phase  function  is 


i.  r  2  2 ' 

+  Lz1  ) 


"  AZJ 


<t>(P1>z1) 


(8.2) 


Let  us  abbreviate  the  derivative  x3/9x^  +  y3/9yi  (x  and  y  are  unit  vectors) 
by  the  symbol  The  stationary-phase  points  for  given  r  and  z\  are  then 
determined  by  VT  ¥  =  0,  which  yields 


P  -  P 


r  -  r. 


7  *  \  \  *<ri> 


(8.3) 


There  are  two  stationary  points,  again,  as  in  Section  6  for  (6.2),  and  the 
backscatter  point  is  ignored  on  the  basis  of  the  discussion  in  Section  5,  which 
can  be  repeated  here  as  an  argument  based  on  terms  (5.1)  except  (i) :  we 
deal  with  6e(p]^,zj)  rather  than  with  6e(Km,zm),  and  (ii):  the  wavenumbers 
kf  and  k^,  are  to  be  replaced  by 


1 

i  — x  .  2 

1/2  ( 

li 

4-i 

1  +  (Ap1/Az1) 

t 

1 

i  — w  ,  2 

1/2  1 

k,  =  -  k  < 
b  1 

j  1  +  (Ap1/Az1) 

+  1i 

(8.4) 


where  (8.3)  is  utilized  for  Ap^/Azp  The  rest  of  the  argument  is  similar. 

We  note  that  the  left-hand  side  of  (8.3)  is  a  direction  sine  independent  of 
location^:  because  the  righthand  side  of  (8.3)  depends  only  on  rp  Hence  we 
name  it  0j .  In  particular,  we  note  chat  (8.3)  holds  in  the  limit  r^  •*  r, 
and  hence  we  obtain  the  ray  equation. 

-v 

"*■  _  dp  1  „  , 

"  d^  =  k  VT  *(r)  ’ 


0 


(8.5) 


where  ds  is  the  length  segment  of  d|r  -  rjJ  in  the  limiting  case  that  -*■  r. 
In  order  to  obtain  (7.2)  from  (8.1)  we  must  be  able  to  set  T  »  $  at  the 
stationary  points.  This  requires  kfL  «  1,  so  we  obtain  the  condition 


kL  <8  >«  1 


(8.6) 


One  should  compare  this  condition  to  (7.3).  Both  conditions  yield  the  Moliere 
approximation,  hence  they  should  amount  to  the  same  thing.  We  will  utilize 
this  to  specify  what  is  meant  by  <(60)2>  in  ^7.3).  The  averaging  in  (8.6)  is 
defined  very  precisely  because  we  know  that  6  is  given  by  (8.5),  and  <f>  is 
given  by  (7.2)  when  (8.6)  holds.  It  is  self-consistent  to  compute <02>  using 
(8.5)  and  (7.2)  and  then  to  apply  it  to  (8.6).  A  simple  calculation [12]  yields 


:02>=  f  dKK3*(K)  »  2L<m1/3LD"2/3 


(8.7) 


Now  compare  (7.3)  to  (8.6)  with  <6  >  given  by  (8.7).  If  these  two  conditions 
are  identical,  then  it  follows  that  o*  <(60)2>  must  be  identified  with  <0^>  . 
However  we  can  replace  a  factor  by  k^  (60)2  in  (8.7).  thus 


o2  <(60) 2>  <02>  = 


oo 

/ 


dKKI-(K)  (60) ‘ 


(8.8) 


2  2 
Note  that  if  we  left  off  (60)  in  the  integral  in  (8.8)  we  would  obtain  a 

by  Itself!  Hence  the  meaning  of  c<2  <  (60 ) 2  >  appears  to  be  that  one  should 

consider  as  an  integral  over  K  (or  over  60  =  K/k)  and  weight  it  by  (69)2 

to  get  the  average.  For  the  Moliere  approximation  we  have 


2  2 


/dKK4> (K)  «  0.4  k2L  Le2 

0 


a2  <602>  =  -  J  dKK$>(K)  (K/k)2  *  2Lk„1/3L„  2/3  e2 


8ir 


m  o 


<(66)2>  =  a2  <662>  /a2  ^  (k  L  )1/3/(kL  )2, 

mo  o 


(8.9) 


and  thus  we  have  defined  the  averaging  procedure  of  the  development  in  Section 
6  regarding  (6.4)  -  (6.6)  quite  precisely. 

Finally  we  can  utilize  (8./)  or  (8.9)  to  yield  (7.3)  or  (8.6)  more  precisely: 
the  limiting  condition  on  the  Moliere  approximation  is  therefore 


21 


(8.10) 


kL2K  1/3L  '2/3e2  «  1 


It  will  be  discussed  in  context  in  Section  19,  where  a  helpful  diagram  will  be 
given  to  illustrate  it  and  similar  conditions. 

As  a  last  comment:  when  we  do  not  drop  the  exponential  function  in  (6.10)  we 
have  to  replace  in  (8.9)  by 

c2  =  f  dKK*(K)  {— } 

o 

where  { — }  is  a  filter  function,  of  magnitude  unity  at  most,  which  filters  out 
some  of  the  integrand.  Hence  (8.9)  is  an  upper  bound  on  our  averaging  proce¬ 
dure;  it  gives  greater  weight  to  large  66  (or  K/k)  than  (8.11)  would.  This 
will  be  an  extremely  useful  property i 
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9.  STATISTICAL  TREATMENT  OF  B:  INTRODUCTION 


So  far,  we  have  found  a  small-angle  scattering  approximation  to  each  term  Bn 

of  the  Bom  series  B  =  1  +  B^  +  B2  +  - .  The  general  term  is  now  given  by 

(6.10)  at  the  observation  point  r  =  (0,L) .  The  expression  for  r  =  (p,L)  is 
easily  observed  from  (3.4b)  -  (3.7)  to  be 


B 

n 


-A  -A  l  £  / 


dz  6e(K  ,z  ) 
m  mm 


-10  (z  -  -  z)  /2k  -  IK.p 
Til  m— 1  m  m 


(9.1) 


A  second  form,  useful  to  us  at  a  later  stage,  is  found  by  replacing  Qm  by 
Km  +  Qnj+1  in  (9.1),  and  rearranging  the  exponential  terms, 
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)/k  -  iK  .p 
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(9.2) 


We  obtained  the  WKB  approximation  by  ignoring  the  first  two  parts  of  the  ex¬ 
ponent  in  (9.2).  The  entire  subject  of  wave  propagation  statistics  in  turbulent 
air  is  determined  by  how  we  will  handle  these  exponential  factors.  It  is  the 
cardinal  problem  in  all  following  development.  Let  us  name 


-iK  2(L  -z  ) /2k  -iK  *Q  (L  -z  )/k 
mm  m  xm+l  ra 

e  e 


.2  ) 
m 


(9.3) 


We  will  eventually  give  the  following  three  developments: 


(a)  In  the  radiowave  regime  we  normally  are  interested  in  propagation 
paths  L  >  kL0^.  In  that  case,  both  factors  of  (9.3)  oscillate,  and 
B  is  the  sum  of  a  constant  and  a  Rayleigh-distributed  field. 

(b)  In  the  optical  regime  we  have  k l0^  <  L  <<  kL0^,  and  the  second  factor 
of  (9.3)  oscillates  less  than  the  first.  Consequently  the  second  can 
be  ignored  when  the  first  is  marginal,  but  when  the  second  factor  is 
dropped  we  can  sum  Bn  of  (9.2)  to  obtain  ch^  Rytov  approximation 
[using  (7.1)  again]. 


B  =  exp 
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/  d2./c 


d  K  6e (K,z)e 


-iK  (L  -  z) /2k  -  iK.p 


(9.4) 
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Act"°i ly,  this  is  not  accurate:  what  we  obtain  above  is  only  a 
factor  c£  proportionality.  The  actual  result  insures  energy  conser¬ 
ve  t_on  whereas  this  factor  does  not. 

( c )  In  the  same  optical  regime  as  under  (b) ,  but  with  important  oscilla¬ 
tions  from  both  factors  of  (9.3)  there  are  to  date  no  analytical 
approximations  that  can  be  trusted  for  B.  Nevertheless,  expressions 
for  <  B>  and  for  <  B(l)  B*(2)>  at  points  $]_  and  1?2  can  be  given, 
and  they  are  identical  to  those  obtained  from  the  WKB  approximation. 
This  regime  is  the  "stong  optics"  regime  with  L  extending  beyond  the 
irradiance- variance  saturation  distance. 

However,  in  order  to  develop  the  theory  —  particularly  for  (a)  —  it  is  essential 
to  discuss  the  treatment  of  A  in  terms  of  statistical  moments.  The  treatment 
of  <  lN>  ,  and  <  B(l)  B*(2)>  differs  slightly  from  that  of  <  B>  where  we  must 
compute  <  Bn> .  The  basic  statistical  problem  lies,  in  the  latter  case,  in 
treating  the  n-point  correlation 


<  6e(l)  F  (1) - 6e(n)  F  (n)> 

n  n 


(0.4) 


where  we  have  used  an  obvious  short-hand  notation  for  the  arguments  of  6c  and 

V 

The  treatment  of  <  1N>  and  <  B(l)  B*(2)>  is  different  only  in  that  we  multiply 
several  terms  uf  B  and/or  B*  together  and  then  take  the  average.  Consider  as 
an  illustration  the  term  <  Bp  Bq*>  of  <  I>  ,  ami  let  n  ■  p  +  q.  We  can  organise 
the  double  multiple  integral 


P-1 


q-1 


j dzx - J  dzp  I"  dz1' - j'  dzq  <  6e(l)  F?(l)  — 6e(q)F(i*(q)> 


(9.5) 


into  n!/p!q!  ordered  integrals 


L  n-1 

/  d‘i'  “  /  K'  <  S'«>  Fn(1>  ” 


-  6e(n)  Fn(n)> 


(9.6) 


where  Z]/ '  —  znM  is  a  permutation  of  z^ - Zp,  z\' - Zq'  that  leaves  the 

order  of  z\  >  Z2  >  -  >  Zp  and  of  zi'  >  Z2*  -  >  zq'  unchanged  but  mixes  z 

and  z'  coordinates.  Therefore  some  of  the  6^(’)Fn(j)  tn  (9.6;  may  be  complex 
conjugates  (although  not  indicated  in  the  notation  here),  but  the  principal 
idea  is  that  we  have  an  ordered  integrand  in  (9.6)  as  in  (9.4)  where  ^ 

0  <  zn  <  -  <  zi  <  L.  Obviously  this  procedure  can  be  extended  for  <  l  >  . 

Consequently,  the  statistical  problem  is  to  reduce  the  n-point  correlation 
<6f(1) - 6^(n)>  to  a  tractable  form.  We  outline  the  procedure  here,  but 
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will  delve  into  dt tail  in  three  subsequent  sections.  Because  Lp  «  L  [L<j  is  a 
measure  of  the  cort  .lation  of  any  twc  of  the  random  variables  6e(i) ] ,  it  fol¬ 
lows  that  if  any  cluster  of  coordinate  t  z\' - Zp'  chosen  arbitrarily  from 

Zi - zr.  (p  ^  n)  is  v*ell  separated  from  all  other  n-p  coordinates,  then  we 

may  write  to  good  approximation. 


<5e(l) - 5e(n)>  ~  <  6e(l') - 6c(p')> 


x  <  6e(p+l') - 6e(n')>  . 


(9.7) 


The  procedure  can  ba  continued  all  the  way  down  to  two-point  correlations. 

Depending  upon  the  way,  the  z\  -  Zq  coordinates  separate  with  respect  to  the 

length  L0,  the  n-point  correlation  will  separate  into  factors  of  lower-order 

correlations .  We  shall  show  two  properties  of  the  n-point  correlation 

<6^(1) - 6^(n)>  with  L  >  z,  > - ->  z  >0: 

i  n 

(i)  The  only  separations  into  clusters  such  as  in  (9.7)  of  importance 

are  those  for  which  6^(1),  -  6e(n)  do  not  change  in  ordering  from 

left  to  iglit.  I.e.,  for  a  4-point  correlation  we  consider 
<6^(1)  6e(2)>  <  6e(3)  6£(4)>  but  not  the  other  two  2-pcint 
correlation  products.  See  Section  11. 

(il)  The  only  separation  of  ordered  clusters  to  be  considered  is 

<6^(1)  fi£(2)> - <  6e(n~l)  fi£(n)>  ,  a  product  of  consecutive 

2-point  correlations,  hence  n  must  be  even.  See  Section  12. 

«JL> 

In  order  to  develop  these  properties,  we  require  a  model  for  Sc(r)  which  will 
be  discussed  next  in  Section  10. 
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10.  AN  EDDY  MODEL  FOR  THE  PERMITTIVITY  DEVIATION 


In  this  section,  a  turbulent-eddy  model  for  6e  is  developed  that  conveys  some 
of  the  intuitive  features  of  fluid  eddies  in  a  turbulent  liquid,  allows  us  to 
further  separate  the  electromagnetic  scattering  part  from  the  statistical  part 
of  the  problem,  yet  is  sufficiently  flexible  to  contain  the  Kolmogorov  spectrum 
as  a  special  case.  The  model  is 


6e(r,t)  * 


f  -  V*> 


(10.1) 


where  Ki  is  a  function  chosen  to  simulate  a  shape  of  an  eddy,  dependent  on  one 
parameter  £.,  and  H.£(t)  is  a  coordinate  that  represents  a  random  trajectory 
through  tne  turbulent  medium.  Wa  also  define  a  density  function  n  (£)  *  N  (i)/V, 
the  ratio  of  the  number  N(£)  of  eddies  defined  by  l  to  the  volume  V  of  the  ran¬ 
dom  medium.  Finally,  we  define  a  Fourier  transform  (<)  of  £t(r)  by 
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(10.2) 


Let  us  first  compute  <  6e(r)>  and  force  it  to  be  zero.  By  inserting  the  in¬ 
verse  of  (10.2)  into  (10.1)  we  find 


<6e(r)> 


(2n) 


-3 
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R„ 


, 3 
d  k 


-iic 
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(r-R,) 


(10.3) 


where  performs  a  random  trajectory  through  space.  If  th^re  are  no  wavenum- 
oers  k,  effectively,  below  a  minimal  value  <min  »  V-1/3  (any  dimension  of  the 
random  medium,  then  it  is  well  known  from  analogous  developments  in  statistical 
mechanics  that 


~i*<'*  R.  o  _ i 

<e  >=8ttV  6^(k), 


(10.4) 


to  good  approximation;  consequently  (10.3)  reduces  to 


6e(r)>  =  y  £  n n  ^0)  =  J2  -^-nf(0) 

R  l 

l  °° 

=  l^n(e)  rip  (0)  =  J  dt  n(£)  n£(0) 


(10.5) 
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The  last  form  represents  a  transition  from  a  discrete  to  a  continuous  descrip¬ 
tion  of  the  scale  sizes  Z.  If  we  use  it,  then  n(£)  has  the  dimension  of  Z~ ^ 
instead  of  Z~^,  In  order  to  obtain  < 6e >  =0,  we  need  only  assume  n(-fc)  = 
n(i)  and  5jt  “  (A/ |&|)  £l&|  (i.e.,  is  negative  for  negative  l) . 

Next,  we  attempt  to  fit  any  spectral  function  $(K)  of  the  permittivity  deviation 
6e.  To  do  sc,  we  compute  <6e(ri)  6e*(?2)>  .  If  we  insert  (10.1)  then  replace 
Kz  by  no  through  the  inverse  of  (10.2),  then  utilize 


> 


» 


we  obtain 


<6e(r1)  6e(r2)> 


/d3K  £  hp(K)|2  n(£) 
8tt  j  l  * 


_V  _X  V 
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(10.6) 


from  which  it  follows  chat 
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J  dZ  n(z)  I n (<) | 2 

—CO 


(10.7) 


Thus,  we  have  one  option  left:  we  can  choose  n(<)  such  that  we  can  solve 
(10.7)  in  n  (Z)  for  given  $(<).  Note  that  we  prefer  the  continuum  notation  in 
(10.7)  for  functions  of  Z. 

An  extremely  practical  choice  appears  to  be  the  Gaussian 
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3/2  3 
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(10.8) 


It  allows  us  to  split  both  Cj,  and  ri£  into  two  factors,  one  dependent  upon  the 
propagation-direction  coordinate,  and  the  other  upon  the  transverse  coordinate. 
Furthermore,  it  enables  us  to  invert  (10.7)  which  first  becomes: 


♦  OO 


d£ 


Z^  n(£) 


2  2 
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(10.9) 


2 

This  is  a  Laplace  transform  with  respect  to  coordinate  Z  .  Let  us  attempt  the 
inverse  for  the  Kolmogorov  spectrum 
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(10.10) 


*0O  -  15.7  L  3  (1  +  A  2)‘11/6  exp(-K2/K  2) 
o  o  m 

2  o 

By  naming  k“  *  p,  and  *  x,  we  obtain 

<2.)5'J  „  [(2*)1'2]  - 

C+i<r 

2^  /  dp 

C-i“ 


and  by  making  some  obvious  transformations,  we  can  modify  the  bracketed  inte¬ 
gral  into  a  factor  containing  (x  -  k^2)  Lq~Z,  and  an  integral,  independent  of 
this  factor,  which  is  essentially  Hankel's  contour  integral  for  the  inverse  cf 
the  gamma  function  [Reference  (13),  Equation  (6.1.4)].  The  result  is 

2  x5/6  -10/3  -*2/2Lo2 

~T^2)  1  6  ’  00.12) 

m 

for  |t  |  >  /2/k^ 

to  good  approximation  (we  have  ignored  a  factor  exp(l/<m2  L02)  in  the  above 
equation)  .  For  |£]  <  /2/ *cm  we  obtain  n(£.)  *  0,  an  artificiality  due  to  the 
choice  of  a  Gaussian  decay  in  (10.10)  and  to  our  choice  of  (10.8).  Note  that 
n(£)  has  no  scalelengths  for  %  <  0.24  i.0,  a  maximum  very  close  to  this  minimal 
value,  namely  at  0.36  iQi  then  decreases  as  2~10/3  until  the  ’-icinity  of 
£  ^  L0  ^2  is  approached,  and  finally  decays  as  a  Gaussian  beyond  this  region. 
The  model  thus  indicates  that  all  essential  eddy  scale  lengths  are  sandwiched 
between  'v*  2.0  and  *,  ^  L0,  quite  in  accord  with  the  physical  model  giving  rise 
to  the  Kolmogorov  spectrum! 
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11.  ELIMINATION  OF  NON-ORDERED  (BACK-SCATTER)  CLUSTERS 


We  now  return  to  the  end  of  Section  9 ,  and  set  out  to  prove  the  first  of  the 
two  assertions  about  the  n-point  correlation  function.  Specifically,  we  will 
consider  a  term  <Bn>  of  <  B>  ,  as  given  by  (9.1).  The  treatment  of  <  I^> 
or  <  B(l)  B*(2)>  is  wholly  analogous  by  virtue  of  (9.5)  and  (9.6).  However, 
in  (9.1),  we  replace  5e(Km,zm)  through  Fourier  transformation 


6e(K  ,z  ) 
m  m 
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-ik  z 
m  m 


(11.1) 


where  Km  =  (K^.k,,,).  Note  that  we  also  use  the  Greek  letter  n  for  the  transform 
in  (10.2),  but  here  there  is  no  sub-index.  It  will  be  clear  from  the  context 
which  transform  is  implied.  Thus  (for  j?  =  0) ,  we  obtain 
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(11.2) 


It  can  be  seen  from^an  argument  similar  to  that  connected  to  (5.1)  that  we  may 
replace  n(Km)  by  nO^.O)  without  appreciable  error.  Consider  the  dzm  integral 
in  (9.2)  for  the  trivial  simplification  (3  =  0,  and  handle  it  as  we  have  handled 
the  first  integral  in  (5.1).  As  illustrated  in  (5.2)  the  variance  of  the 
integral  is  proportional  to  <t[(Km-  +  kf2)l/2]>  where  in  this  case 
kf  =  Km(Qx,i  +  Qm+i)/2k.  The  point  is  that  (Qm  +  Qm+'^/Zk  Qm/k  1  and  there¬ 
fore  f  <<  Km-  If  we  ignore  kf  in^Q,  we  replace  $[Km,kf)  by  $(Km,0).  Because 
$(km»kf)  is  proportional  to  <|  n(Km,kf)|2>,  and  because  the  very  same  argument 
must  hold  in  (11 .jp  for  the  combined  dzm  dkm  integral,  we  note  that  o(Km,km)  can 
be  replaced  by  o(Km,0)  to  the  same  order  of  approximation.  The  error  is  ob¬ 
viously  0(kf2/Km2)  =  o  (Qm^/k2)  in  each  m-th  factor,  and  is  therefore  cumula¬ 
tively  effective  as  described  in  Section  6. 

The  next  step  is  to  perform  th^  dz:a  integrals  in  (11.2).  To  free  ourselves 
from  the  difficult  boundaries,  we  introduce  an  extra  contour  integral  just  on 
the  positive  side  of  the  imaginary  axis  of  a  new  complex  variable  p.  That  is 
to  say,  we  utilize  the  equalit>, 
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The  conditions  on  the  second  form  in  (11.3)  are  explained  by  the  fact  that 

L  -  E  Az  *  z  .  Hence,  z  >  0  and  Az  >0  is  identical  to  0  <  z  <  -  <z.  <  L. 

m  m  n  n  m  n  1 

The  third  form  rids  us  of  the  condition  L  -  E  Az^j  i  0  and  therefore  allows  us 

to  boi  •  d  each  dAzm  integral  between  0  and  The  reason  is  that  the  contour 

integi.^  (c  is  a  vanishingly  small  positive  increment)  of  (2rrip)~l  exp[p(L  ~  E  Az  )  ] 

is  unity  for  L  -  E  Azm  i  o„  and  zero  otherwise.  Because  the  dAzm  integrations 

are  now  independent,  they  can  be  performed  to  yield, 
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Now  we  introduce  the  eddy  mgdel  of  Section  10.  However,  we  have  to  do  it  in 
(11.2),  where  we  replace  n(k  )  by 
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and  then  utilize  the  last  form  in  (11.4).  Thus,  invering  the  order  cf  integra¬ 
tion  we  obtain 
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where  it  should  be  noted  that  pg,  zg  are  functions  of  index  m  (although  not  so 
labeled) ,  and  that  the  d2Rm  integrals  should  be  to  the  left  of  all  the  factors 

containing  Qx,  Q2 - ,  Qn.  Let  us  rename  *  (pm' ,  z^' ) ,  so  that  primed 

coordinates  represent  the  randomly  moving  centers  of  eddies.  Then,  we  trans¬ 
form  all  kro  to  qm  (Jacobian  is  unity)  which  transforms  -ikm(L  -  %')  to 
-iqm  Azm' ,  where  Azm' ,  =  Zm'  -  Zj„.  The  product  of  the  dkm  integrals  in  (11.5) 
thus  becomes 

30  -iq  Az  ' 

n  -  na  m 

n  /  da  - - - 5 -  =  (11.6) 

m=l  «/  p-iq^  +  1  Qm  /2k 


Note  that  p  has  a  small  positive  real  increment  c.  Consequently  the  denomina¬ 
tors  have  zeroes  just  below  the  real  axis  in  die  complex  qm  planes.  It  there¬ 
fore  follows  that  (11.6)  is  zero  when  any  of  the  Azm'  <  0  because  the  qm  con¬ 
tour  can  then  be  closed  in  the  half  plane  without  poles.  As  a  result  Bn  is 

non-zero  only  when  all  A >  0.  Thus  we  must  have  0  <  zn'  <  -  <  zx*  <  L, 

in  which  case  (11.6)  reduces  to. 
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Upon  inserting  (11.7)  for  the  factors  (11.6)  into  (11.5),  we  can  perform  the 
dp  contour  integral  to  eliminate  the  p  dependence  again,  and  we  obtain 
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This  is  the  result  required  to  prove  the  property  (i)  at  the  end  of  Section  9. 
We  note  that  the  n-point  correlation  in  <  Bn> ,  in  terms  of  (11.8),  is  deter¬ 
mined  by 


<  E 


lie,  *R 

X 

e 


fc.i 


i< 


e 


n 


> 


(11.9) 


where  xm  =  (K^,  (Qm  -  Qm+l)/2k).  As  discussed  in  Section  10  the  Ri,m  'ol- 
low  random  trajectories;  therefore  (11.9)  breaks  up  into  products  of  lower 
m-p<  int  correlations  (m  ^  n)  or  clusters  so  that  any  two  exponentials  in  (11.9) 
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with  %  ^  «  R£tj  belong  in  two  different  clusters  and  any  two  exgonentials  with 
identical.  Rj^i  are  in  the  same  cluster.  However  we  cannot  have  R£>nrf2  = 

and  ^d,m+l  ^  %  m  because  of  the  ordering  condition  0  <  Az^,  <  l  of  (11.8). 
Therefore  the  only  way  that  (11.9)  factors  is  in  the  fashion  described  by  (i) 
at  the  end  of  Section  9,  namely  in  ordered  clusters  so  that  aftej  making  the 
division  into  products  of  separately  bracketed  factors,  all  the  R£jta  are  still 
in  the  same  order  as  in  (11.9). 

In  thus  eliminating  the  n2n-ordered  cluster,  we  have  made  only  one  further 
approximation,  namely  n£(Km,  km)  nj,  (Km,  0).  We  have  proven  that  this  gives 
an  error  of  0  (Q^/k^)  which  has  the  cumulative  effect  discussed  in  Section  6, 
and  overspecified  in  Section  8. 
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12.  REDUCTION  TO  THE  MAIN  BINARY-CORRELATION  PRODUCT 


So  far  we  have  demonstrated  that  the  n-point  correlation  (9.4)  or,  equivalently, 

(11.9)  may  be  written  as  a  sum  of  terms  obtained  by  subdividing  R£  =  R^,  R2, - , 

Rn  into  all  possible  ordered  groups  of  singles,  pairs,  triplets,  etc.  which 

retain  the  sequence  1,  - ,  n,  and  with  averaging  brackets  around  each 

m-tuplet  (which  is  therefore  an  m-point  correlation  or  an  m-th  order  cluster). 
Now  we  will  demonstrate  that  ouly  the  term  consisting  of  n/2  binary  clusters 
is  of  importance. 

Let  us  take  Rj ,  - ,Rjj  in  (11.9)  -  note  that  we  have  dropped  the  extra  index  & 

for  convenience  -  and  subdivide  it  into: 

m^  groups  of  uncorrelated  singles  R^,  ___ 
m2  groups  of  uncorrelated  pairs  Ri  =  R j ,  etc. 

-X.  -X  -X 

Thus,  a  group  with  Rp^  =  Rp2  - - =  yields  a  factor 

<exp[i(ic  ^  + - +  Kpj)*®pj?'>  an^  there  are  mj  factors  which  have,  exactly  j 

wavevectofs  k  in  the  exponential.  Note  that  we  allow  m^  0  at  this  time 
even  though  we  already  know  that  <6e(K,z)>  =  0.  The  condition  on  mj  is 

n 

E  jmj  =  n  (12.1) 

j=l 


The  above  factor  is  thus  a  particular  one  in  a  term  of  (11.5).  There  are  many 
ways  in  which  (11.9)  can  be  split  into  m^  groups  of  1,  m2  groups  of  2,  etc. 
products  of  factors;  the  combinatorial  part  of  the  problem  will  be  dealt  with 
later.  Let  us  reconsider  the  j-th  product  and  rewrite  it  in  simplifies  nota¬ 
tion  by  dropping  the  index  p.  We  invoke  the  two-dimensional  analog  of  (10.4) 
and  find 

1 

<  exp[i  E  V\l> 

m=l 

where  d^  is  the  width  of  the_jnedium  in  a  direction  perpendicular  to  z  (we  let 
drp  -v  «  ultimately),  and  6 2 (^Kpi)  is  the  two-dimensional  Dirac  delta  function. 

Let  us^see  what  +  -  +  •=  0  means.  Consider  (11.8)  and  identify  k^ 

with  (Kflj+i » (Qj^  “  Qm+2^^  '  means  first  of  all  that 

K.  2  =  ^m+2  = - =  K  (lnc,ex  m  r‘ee^  not  be  specified  here:  it  depends  on 

where  this  particular  j-cluster  is  sandwiched  in  (11.9).  It  also  means  that 

Qm+1  =  QnH-2  = - =  Qnrt-j+l>  so  that  =  (Kj»°>  for  10+1  ^  P  k  •  So  if  we 

^elect  from^(11.8)  only  those  factors  and  integrations  pertaining  to 

K,^ » - ,  K^j  (and  dropping  index  m  for  convenience  we  have 
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Note:  It  is  crucial  to  note  that  the  disappearance  of  all  Qm  vectors  from 
(12.3)  is  not  general.  It  holds  for  <Bn>  ,  and  for  <  Bp  Bq*  >  with  p+q  =  n, 
but  it  does  not  hold  for  higher-order  field  correlations.  In  those  other 
cases,  we  obtain  terms  -  ic^  as  well  as  +  Kp  (1  ^  m,p  j)  in  (12.2),  and  as  a 
result  the  wavenumbers  (Q^-i  -  Q^2^2k  r'ot  allu.cancel •  ^Therefore,  we 

really  have  to  add  a  function  of  the  form  exp[i  f(Qi,  - ,  CM  ]  as  a  factor 

of  the  integrand  in  (12.3).  This  is  an  oscillating  factor  with  unit  amplitude, 
hence  we  overestimate  the  integral  if  we  exclude  it.  The  error  thus  made  will 
not  harm  the  results  if  we  utilize  the  overestimate  oi  ly  to  estimate  the 
magnitude  of  terms  ultimately  discarded  in  comparison  to  terms  kept  (where 
this  error  need  not  be  made) . 


In  order  to  continue  the  estimate  (exact,  rather  than  an  estimate,  only  for 
<B>_iand  <B(1)B*(2)>  ),  we  utilize  (10.8),  and  the  degeneracy  of  the  sum 
over  R j ,  i . e . 


£  f(i)  =  LdT2  f  dl  n 


U)f  (£). 


(12.4) 
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This  property  is  the  same  as  the  one  used  in  (10.5)  with  V  =  Ldf  .  Hence, 
(12.3)  becomes 
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We  free  the  d  integration  from  interdependence  writing  the  delta  func¬ 
tion  as  an  integral  of  (4tt2)”1  exp[i(Kj^  +  —  +  Kj)*p]  over  d^p .  Then,  we 
utilize 
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to  obtain 
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Because  there  are  mj  factors  of  (11.8)  pertaining  to  a  j  cluster,  and  because 
we  allow  1  ■£  j  4  n,  we  obtain  for  one  factorization  characterized  by 
{mj}  =  {mi,  m2  - ,  1%} 


j+2 

d£  n(A)£ 


(12.8) 


This  is,  we  reiterate,  correct  for  <B>  and  tor.  <B(1)B*(2)>  ,  but  an  over¬ 
estimate  for  other  statistics  of  B.  However  it  will  serve  as  a  basis  for 
estimating  the  effect  of  decomposing  the  n-point  correlation  into  the 
ml,  m2»  - *  ®n  product  of  mj  factors  of  j  clusters.  Let  us  denote  this  de¬ 

composition  by  {mj}  as  indicated  in  (12. 8> . 

First  of  all  we  note  that  {mj}  cannot  contain  any  mj  t  0  for  odd  j.  The 
reason  is  that  n(-2.)  =  n(£) ,  hence  (12.8)  is  zero  for  any  occurrence  of  an 
odd  j.  Therefore  only  even  clusters  ccntribute. 

The  rest  of  this  section  will  be  devoted  to  demonstrating  the  predominance  of 
the  binary-cluster  decomposition.  Because  n  must  be  even,  we  henceforth  set 
n  =  2M.  The  binary  cluster  is  characterized  by  m2  =  M,  all  other  mj  =  0. 

Let  us  now  estimate  (12.8).  By  utilizing  (10.12),  or  other  physical  spectra, 
we  note  that 


J  d£  n(£)£j+2  «  L/  1 


Consequently,  we  note  that  (12.8)  is  estimated  by 
2  2  M 

(-k  LL  e  )  (L  /L)  with  m  =  >, 

0  0  J 

1=1 


(12.9) 


(12.10) 


We  have  made  use  of  (.12.1)  in  arriving  at  this  estimate.  Because  L0  <<  L,  it 
follows  that  (12.10)  is  maximal  for  m  =  £  nij  maximal.  However,  m  is  subject 
to  the  restriction  (12.1),  and  therefore  m  is  maximal  for  m2  =  M,  all  other 
mj  =  0,  i.e.,  for  m  =  M.  We  have  proven  that  the  binary-cluster  decomposition 
is  dominant,  but  not  that  the  cumulative  effect  of  other  decompositions  can  be 
neglected.  This  is  easily  done  with  (12.10): 

Note  that  m  =  E  mj  is  the  number  of  delta  functions  occurring  in  (12.3)  factors 
of  (11.8)  for  n  =  2M.  Allowing  only  even  j,  there  are  (M-l) ! /m! (M-m-1) ! 
different  ways  of  partitioning  the  2M  coordinates  into  m  consecutive  groups  of 
even  coordinates  or,  equivalently,  into  m  delta  functions.  The  error  estimate 
based  on  (12.10)  is 
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(12.11) 
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Consequently,  the  error  in  discarding  other  than  the  binary-cluster  decomposi¬ 
tion  is 


0(MLq/L)  =  0{a\QIV) 


(12.12) 


2 

where  a  is  the  expansion  parameter  discussed  in  Sections  4,  6,  and  8.  Note 
that  (12.12)  is  identical  to  the  first  of  the  errors  (6.6),  hence  the  present 
development  is  also  important  in  understanding  the  nature  of  the  approximations 
of  Section  6,  in  particular  the  discussion  of  (6.5).  Note  also  the  relation¬ 
ship  of  (12.10)  to  (8.9). 

2 

A  special  case  of  (12.12)  follows  for  <  B  >  ,  because  a  can  be  identified 
with  k2LLce2,  see  (12.10).  Thus  the  predominance  of  the  binary-cluster  term 
and,  concomitantly,  small-angle  scattering  must  require  <j 2 L0/L  <<  1,  or 
k2L02e2  <<  1.  Later,  we  gather  all  the  approximations  together,  but  here  we 
note  that  we  name  this  as  the  condition  for  the  binary-cluster  expansion.  In 
othj^rjtfords,  the  Rendition  k2LQ2e2  <;-  1  is  sufficient  for  decomposing 
<6e(K1,z1) - 6e(K2m,Z2M^>  as  <  :(K1,z1)6e  (K2,z2)  > - 

<  <5e(K2M-l»z2M-l)  SeOC^.z^j^  in  jrroducts  of  terms  of  the  Born  series, 
provided  the  other  small-angle  scat;ering  approximations  are  met. 
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13.  STATISTICAL  TREATMENT  OF  B:  DIAGRAMS 


After  the  excursion  of  Sections  10,  11,  and  12,  we  are  ready  to  return  to 
Section  9,  and  resume  the  further  statistical  treatment  of  B.  Let  us  summarize 
the  status  at  this  point: 
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Furthermore,  we  may  replace  the  2M-point  correlation  in  any  statistic  of  B  by 
the  ordered  binary-correlation  product,  i.e., 


<6e(l) 


6e(2M)> 


2M-1 

j1  <6e(m)6e(nri-l)> 
m=l 


(13.2) 


The  remaining  problem  -  aside  from  actual  calculation  of  terts  -  is  one  of 
bookkeeping.  It  is  easy,  of  course,  for  <B>,  but  in  computing  <B(1)B*(2)>, 
it  is  already  more  complicated  because  some  of  the  factors  <5e(m)  in  B(l) 
combine  with  other  contiguous  factors  in  B(l),  but  other  ones  combine  with 
factors  6e*(n)  of  B*(2).  The  complications  grov?  rapidly  when  computing  higher- 
order  moments  of  I,  <  1^  > . 


In  order  that  the  bookkeeping  problem  becomes  tractable,  we  introduce  diagrams. 
As  an  example  consider  a  term  of  <  I^>  =  IQ^  <(BB*)^>  of  order  2M  in  6e. 

I.e,  there  are  N  factors  Bn^,  and  N  factors  B*mj  (where  ni  and  mj  are  equiva¬ 
lent  notations  for  n^  and  m j ,  respectively),  such  that 


N 

E  K  +  m.)  =  2M  (13.3) 

i=l  1 


We  will  use  a  diagram  to  depict  the  term  <  BniB£i  -  BmNBmN'>  *  Actually, 

there  will  be  many  diagrams  for  one  term  characterized  by  the  unique  set  of 

indices  {n-^,mj},  i.e.,  for  a  unique  choice  of  n^,n2, - n^,  m^.m^, - mN 

such  that  (13.3)  is  satisfied.  However,  this  will  turn  out  to  be  an  aid 
rather  than  a  liability!  The  rules  for  a  diagram  are  very  simple: 

(i)  Each  factor  Bn^  is  represented  by  a  full  horizontal  line  (between 
z=0  and  z=L) ;  an  axis. 

(ii)  Each  factor  Bmj  is  represented  by  a  dashed  horizontal  line;  a 
conjugate  axis. 
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(iii)  Each  <Se(m)5&(nrt-l)>  or  <  s£*(m)  6e*(m+l)  >  ,  corresponding  to 

a  correlation  in  one  or  in  one  factor,  is  represented  by  a 
snail  circle  (or  a  large  dot),  which  we  shall  call  a  "bead",  on  the 
corresponding  horizontal  line. 

(iv)  Each  <  62(a)  S£(nri-1)>  or  <6e  (m)62*(mfl)>  or  <  6e(tn)62*(nri-l)  > 
or  <  62*(n)62(ort-l)  >  corresponding  to  correlations  between  two 
different  horizontal  lines  will  be  represented  by  a  vertical 
line  connecting  them;  we  name  the  connection  a  "rung". 

We  wish  to  show  that,  conversely,  any  diagram  of  the  above  type  containing 
beads  and  rungs  corresponds  uniquely  to  a  contribution  to  <I">  (for  this 
example) .  To  do  so ,  we  compute  the  contribution  of  typi cal  beads  and  rungs : 

A  Bead  Contribution:  The  contribution  of  <6e(m)62(m+l) >  as  in  (iii)  of  the 
above  diagram  rules  is  given  by 


z  ,  z 

m-1  m 
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(13.4) 
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but  we  must  keep  in  mind  that  both  d  Kjn  and  d  integrals  are  also  over 

preceding  factors,  possibly .Jbecause  these  may  depend  on  Q^,  Qm_i,  etc.,  all 
of  wh^ch  c^n  contain  K_  and  K^jJsee  (13 . 1)  j .  By  utilizing  (5.2),  we  note 
that  Km  +  =  0.  That  implies^that  (13.4)  is  ngw  independent  of  all 

preceding  factors  because  K^j  and  drop  out  of  Qm,Qm_i,^etc.  By  further 

use  of  (5.2),  and  by  noting  that  K,^  Aznri.j_/2k  «  1  and  Km’Qnri-2  ^znr+d/k  <<:  1> 
we  obtain 


where  is  the  integral  scale,  defined  implicitly  in  (13.5)  by  the  second  form. 
The  definition  (see  Reference  [12],  for  example)  ’  ields  =  0.188Lo.  Note  the 
correspondence  with  a  single  factor  of  (12.10);  we  have,  of  course,  done  the 
same  calculation  there  by  another  method.  We  should  note  that  we  have  as 
integrand  in  the  integial  (13.5)  possibly  a  function  of  zm  because 
^2^m»zm_zm+l)  binds  the  dzmf2_^nte8rab  upper  bound  to  the  immediate  vicinity 
of  zm  (because  z^f^  «  zm) ,  and  therefore  the  result  of  the  factors  beyond 


may  be  a  function  of  zm.  Consequently  the  dzm  integral  in  the  last  of 
(13.5)  is  left  open  ended. 

A  Rung  Contribution  Between  Bn  and  %*.-  In  this  case  we  consider  two  consecu¬ 
tive  6e  fac-ors  [consecutive  in  the  sense  discussed  in  connection  with  (9.5) 
and  (9.6)]  belonging  to  a  Bq  and  a  Bra*  axis  respectively.  Their  contribution 
is. 
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where  the  primed  coordinates  are  associated  with  the  Bjj*  factor.  Note  that 
Zp  and  Zq  are  consecutive  z  coordinates.  We  utilize  (5_^2)^for  the  above  two- 
point  correlation,  noting  that  it  contains  a  factor  A2(Kp~Kq)  so  that  Kq=Kp 
in  (13.6).  After  some  algebraic  manipulation,  and  judicious  use  of  tlje  re¬ 
striction  of  Zp-Zq  to  distances  at  best  several  times  L0  [because  4>2(Kp>zp"zq) 
falls  off  rapidly  otherwise],  (13.6)  reduces  to: 


k2e2 
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where  Ap  =  p-p  .  We  have  assumed  that  there  is  no  zq_i  coordinate  between 
zp-l  and  zp»  otherwise  the  dZp  integration  would  be  Bounded  by  Zq_^.  It  is 
difficult  to  maintain  a  consistent  notation  because  it  would  become  unwieldy 
rapidly.  For  instance,  Qp+i  is  the  sum  of  n-p  wavevectors  of  Bn,  but  *s 

the  sum  of  m-q  wavevectors  of  B™*.  However,  barring  such  inherent  implications 


of  the  notation,  the  result  (13.7)  can  be  read  off  from  the  diagram  if  we  draw 
the  rung  at  location  z=Zp  between  Bn  and  Bm*.  The  wavenumber  Qp+i  is  the  sum 
of  all  rungs  with  one  end  on  Bn,  and  Qq+i  of  all  rungs  with  one  end  on  Bm 
both  groups  to  the  left  of  the  rung  at  z=Zp.  The  integrand  of  (13.7)  may  con¬ 
tain  other  factors  from  rungs  to  the  night  of  the  z  ^  rung  because  these  may 
be  functions  of  Kp*  and  also  from  rungs  to  the  left  Decause  these  may  contribute 
a  function  of  Zp  (the  upper  bound  of  the  dzp+i  integral  is  Zp) .  So  care  must 
be  taken  not  to  consider  contribution  (13.7)  in  isolation  of  o^her  factors. 

Note  that^there  is  no  reason  to  prefer  the  ordering^such  that  Qp+i  -  Qq+i  in- 
stead  of  0q+^  -  Qp+i  occurs  in  (13.7).  Because  K^Kp,  the  other  ordering 
obviously  changes  the  sign  in  both  exponentials.  ‘However,  we  adhere  to  the 
convention  that  the  Q  and  the  p  of  the  top  line  is  to  come  first  (minus  signs 
in  the  exponent  for  top  Bn,  plus  signs  for  top  Bm*) .  Symmetry  properties  will 
eliminate  all  ambiguities. 
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A  Rung  Contribution  Between  Bn  and  (or  Bn*  and  B^):  This  case  is  treated 
nearly  identically  to  the  previous^one,  except  for  a  factor  i2  =  -l  extra  in 
(13. 6>,  replacement  of  -K^p’  by  +K^*p",  and  no  complex  conjugate  factors.  In 
this  case  Kq=K_,  and  if  we  denote  the  coordinates  of  the  lower  B™  by  double 
primes ,  we  then  find. 


k2e2  ^ 
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for  thg  contribution  oi  a  Bj^  rung,  and  the  complex  conjugate,  of  course, 
for  Bq  Bq  rungs.  Note  that  there  is  an  extra  factor,  compared  to  (13.7). 
It  arises  from  the  first  exponential  of  F  given  ir  (13.1).  In  the  cases  of 
(13.7)  and  (13.5)  this  exponential,  cancels  out  to  good  approximation. 


The  rules  (13.5),  (13.7),  and  (13.8)  constitute  a  unique  interpretation  of 
any  diagram.  A  bead  or  rung  at  zp  (0  <  z  <  L)  obviously  defines  all 
Zp ‘dependent  factors  and  integrals  in  these  three  diagrams.  The  Q  wave- 
vectors  can  be  read  from  all  features  to  the  left  that  are  connected  to  either 
or  both  of  the  axes  defined  by  the  zp  feature.  Finally,  Zp_i  is  the  integra¬ 
tion  variable  for  the  next  feature  to  the  right.  We  give  some  examples  in 
Figure  2  f or  < I  >  .  For  Figure  2a  we  note  that  (13.7)  holds  for  the  feature 
31-  zp  v£satures  to  the  right  have  nit  Jr^en  sketched  in^.  The^Q  vectors  are: 

|p+3  +  I?*4  anlQq+l  *  K+2  '*  Xp+3-  N°te  that  KL.3  *  L+3  subtracts 
out  in  Qp+i  -  Q|+1*  The  reason  is,  ofPcourse,  tha^  the.  p+3  feature  is  be¬ 
tween  the  same  two  rungs.  It  is  easily  seen  that  Qp+1  -  Q'+1  does  not  contain 
K-vectors  of  rungs  between  the  same  axes  as  the  Zp  rung  when  two  mutually 
conjugate  axes  are  connected,  furthermore,  beads  do  not  contribute  to  the 
Q  vectors  of  rungs.  Figure  2b  is  nearly  the  same  as  Figure  2a,  except  that 
Jjhe  zp  feature^is  a  rung^of  the^type  described  b£  (13.8)^  Consequently , 

nP+1  *  IR*3  +  ^p+4’  an^  Q?xl  =  KP+4-  Note  that  Kp+4  =  “  xf+4»  hence 
jP+1  ~  ^q+l  =  Kp+3  +  2Kp+4»  Unfortunately,  other  rungs  bitween  the  two  axes 
do  not  cancel  out  their  K-vectors  in  this  case.  The  other  parts  of  the  in¬ 
tegrals  are  easily  and  trivially  written  down. 


tarms  of  °<Lo/D  have  been  discarded  in  the  exponential  factors  of 
(13.7)  and  (13.8),  as  well  as  in  (13.5)  where  no  exponential  factor  is  left. 
This  type  of  error  has  been  discussed  in  Section  4,  and  cumulative  error  is 
the  same  as  in  (12.12). 


Finally,  we  return  to  the  bookkeeping  problem.  The  results  indicate  that  the 
diagrams  represent  a  much  easier  bookkeeping  than  the  products  of  terms  in 
the  Born  series.  ^To  clarify  this,  consider  contributions  to  <B(1)B*V2)> 
proportional  to  €  .  These  come  from  <  B9(1)B  *(2)>.  <  Bi  (l'iBo* > 

<83(1)6!  (2)>  ,  <B4(1)>  ,  and  <TB4*(2^>  .  On  the  otheAandAey  are  given 
by  the  sum  of  all  diagrams  with  exactly  two  features,  to  wit  2  beads,  2  rungs 
ana  1  bead  pius  1  rung,  in  all  possible  sequences.  Not  only  is  this  pictorially 
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Two  examples  of 


2 

<X  > 


diagrams . 


simple;  it  will  also  yield  the  usual  graph-summation  rules  that  make  it 
possible  to  sum  up  series  of  terms  that  appear  to  have  a  complicated  book¬ 
keeping  otherwise.  Quite  specifically,  for  any  statistic  of  B,  the  sum  of 
all  products  of  Born-series  terms  can  be  represented  as  above  by  a  sum  over 
all  topologically  different  diagrams! 


14.  SUMMATION  OF  BEAD  COMPONENTS:  COHERENT  FIELD  <  B> 


One  of  the  first  selective-summation  advantages  of  diagram  techniques  is  the 
elimination  of  all  coherent-wave  effects  in  statistics  -of  B.  That  is  to  say, 
all  bead  contributions  can  be  summed  out,  which  means  that  we  need  only  con¬ 
sider  rung  diagrams  afterwards.  No  correlations  inside  one  factor  Bn  or  Bn,* 
are  left.  We  spoke  of  "coherent-wave"  effects  for  such  correlations  because 
their  sum  ;  mounts  to  the  coherent  or  average  field  <  B>  . 

It  is  very  easy  to  sum  out  all  beads.  Consider  one  bead  at  %  as  in  (13.5), 
and  a  rung  at  z^-i  and  at  z^i,  no  matter  what  statistic  of  B  we  consider 
(except  <  B>).  In  general  we  will  have 

z  „  z  ,  z 

m-2  m-1  m 

/  dVl  £  <Vl>  (-k2<20  /  dzm  /  dVl  8  <Vl>  <U-1> 

o  oo 


If  we  keep  features  m-1  and  m+1  fixed,  then  the  bead  contribution  to  (14.1) 
[where  we  use  f  (%_].)  and  g  (zm_i)  simply  to  denote  contributions  such  as 
(13.7)  and  (13.8)  for  features  (m-i)  and  (nri-1)]  can  be  written, 

Vl 

-k2'2“i  /  dz*  -  -k2'2,i  <Vi  -  Vi> 

Zm+1 


If  there  were  no  bead  at  all, 
upper  bound  of  dz^^  would  be 

have  p  beads  at  zj',  - ,  Zp' 

instead  of  (14.2), 


then  (14,2)  would  be  replaced 
z^l  in  (14.1).  It  is  easily 
with  Zm+1  <  zp'  <  -  <  zi'  < 


by  unity,  and  the 
seen  that  if  we 
zm-l»  we  obtain 


(-  k2<\) 


Zm+1 


(14.3) 


To  add  up  the  effect  of  any  number  of  beads  at  zm,  we  thus  sum  (14.3)  over  all 
p  from  p  =  0  to  p  =  00  to  obtain 


exp 


2  2 

k  e  i. 


(14.4) 
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where  Az  is  the  interval  between  features  m-1  and  mfl.  Thus  if  we  draw  a 
diagram  with  rungs,  we  can  generate  all  diagrams  w’th  rungs  in  the  same  loca¬ 
tion,  but  any  number  of  beads,  by  inserting  factors  (14.4)  for  all  Az,  includ¬ 
ing  the  cases  where  z^-i  =  I.  and  Zm+i  =0.  We  can  multiply  all  the  factors 
(14.4)  for  one  B  or  B*  factor  together  and  we  note  that  each  B  and  each  B*  fac¬ 
tor  in  any  statistic  of  B  (say  <I^> )  contains  a  factor 

exp  ^-k2e2£i  L^,  (14.5) 


and  the  resulting  factors  are  describ°d  by  the  sum  of  all  diagrams  with  rungs 
contributing  factors  such  as  (13.7)  and  (13.8).  This  represents  a  simplifica¬ 
tion  since  it  is  much  easier  to  keep  track  of  only  one  type  of  topological 
feature.  In  particular  (14.5)  is  the  entire  answer  for  <B>  because  there 
are  no  rungs  at  all.  Thus,  a  corollary  of  (14.5)  is 


<B>  =  exp  (-  a  L) 

_  ,  2  2. 
a  =  k  e 

i 


(14.61 


Note  that  we  can  also  obtain  (14.6)  from  the  average  of  (7.2)  by  utilizing  the 
development  leading  to  (7.4).  However,  we  have  already  noted  that  (7.2)  is  of 
restricted  validity  in  terms  of  (7.5),  where  cj2  =  a.  The  above  result  (14. 
is  derived  under  much  less  restrictive  assumptions:  only  the  cluster  assumption 
(12.12),  the  small-angle  assumption  (6.9),  and  the  sagittal  approximation  (6.11) 
have  been  assumed.  Let  us  examine  these  with  the  method  described  at  the  er.d 
of  Section  8,  i.e.,  by  interpreting  60  as  K/k-,  and  inserting  that  into 
o2  =  a  =  k2e  as  given  by  the  second  form  in  (13.5): 


(i)  The  binary -cluster  approximation:  As  discussed  in  (12.12)  we  mean 
the  approximation  obtained  by  ignoring  errors  of  0  (A^,c)  where 
Abe  %  °2  Lo/L.  In  this  case  we  simply  set  a2  ^  k2e2LDL  to  obtain 
Abe  ^  k2L2e2. 


(ii)  The  small-angle  scattering  approximation:  This  arises  by  ignoring 

terms  of  0  (Qm2/k2) .  I.e.,  the  error  is  Aq  ^  <Qm2/k2>  'v  q2  <(60) 2> 

as  discussed  in  Section  u.  Following  the  just  given  recipe,  we 
obt  lin. 


.  Le 

A  - 

0  16tt 


/ 


dK  K3<i>(K)  *  2L  k  1/3  C  2 
m  n 


L  K  1/3  L  -2/V 

m  o 


(14.7) 
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(iii)  The  sagittal  approximation:  Ignoring  terms  of  order  given  by  (6.11) 
requires  us  to  be  more  precise*  Strictly  speaking,  the  expansion 
(6.11)  comes^from^OL  -  Qfji/k2)l/2  -  (1  -  Qmfl/k2)l/2j  an(j  by  substi¬ 
tuting  Qm  “  Km  +  Qmfi,  we  note  that  it  arises  from 

[o  -  -  (4 + 2K  •  w/k2]1/2  -C1  -  '^i7"2)172 

The  noteworthy  feature  of  this  is  that  the  0  (K^/k^)  term  is  (K^  +  21^  • 
provided  Ag  «  1.  Therefore,  we  can  estimate  the  cumulative  sagittal  error  by 
considering  the  effect  of  an  exponential  term  L  Qmfl/k^,  i.e.,  by  estimating 
when 


/ 


dK  K$(K) 


o 


2  2,3  ® 

iLK  <  Q  >  /k  r 

e  *  J  dK  K4>(K) 


We  utilize  the  spectral  form  $(K)  *  (1  +  K2L02)  to  set  x  =  K2L02,  hence 

to  see  when 


00 

f  dx  (l  +  x2) 


2  3  2  00 

-11/6  ixL  <Q  >/k  Lq 

e  ^ 


*  / d*  ( 1  +  *2) 


-11/6 


2  3  2 

rbviously,  th's  is  the  case  when  As  <<  1  if  we  define  As  'v  L  <Q  > /k  Lq  'v 
^  L  A0/kLo2,  thus 


.  .2  1/3  T  -8/3  ,-12 

A  'v  L  k  L  k  e 
s  mo 


(14.8) 


Thus,  there  are  three  important  cumulative 
and  ns.  We  return  to  these  in  Section  19. 
may  be  overestimated  with  respect  to  their 
was  discussed  in  connection  with  (8.11). 


rs  to  keep  track  of:  Avc,  Ag, 
ihould  be  noted  that  the  errors 
.uulative  effect  in  the  sense  that 
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15.  THE  MUTUAL  COHERENCE  FACTOR  <  B(1)B*(2)> 


The  methods  of  Sections  (13)  and  (14)  are  easily  extended  for  calculation  of 
the  mutual  coherence  factor  <B(1)B*(2)>  where  (1)  stands  for  coordinate 

ri  =  (yp,  L),  and  (2)  for  r2  =  (-  j  p ,  I). 

From  the  preceding  it  follows  that  <B(1)B*(2)>  is  given  by  exp(-2aL)  times 
the  sum  of  all  diagrams  with  2  axes  and  m  rungs  (0  4  m) .  The  only  type  of 
rung  that  can  occur  ijj.  that  with  contribution  (13.7) .  £t  is  immediately 
apparent  that  Qp+i  -  Qq+i  =  0  because  q=p  and  all  =  Kp.  The  contribution 
of  the  p-th  rung  is  thus 


dKK$ (K)  Jq (Kp )  =  2a (p) 


o 


(15.1) 


and  hence  we  find  that  the  m-rung  diagram  contributes  in  total 
e“2aL  [2a(p)L]m/m! 


(15.2) 


where  a(p)  is  defined  in  (15.1).  Obviously  a(p)  +  a  as  p  +  0.  Upon  summing 
over  all  m  0  we  find, 

<  B(1)B*(2)  >  =  exp[-2Ao.(p)L] 

Aa(p)  =  a(0)  -a(p)  (15.3) 

00 

kV  f 

=  fjrp  J  dKX*(K)  (1  -  J0(Kp)] 

o 

The  result  is  quite  well  known  and  we  will  not  discuss  it  further,  except  to 

point  out  that  it  is  subject  to  the  same  constraints  &bc  <<  1  and  Aq  <<  1  as 

(14.6).  A  corollary  of  (15.3)  is  energy  conservation:  for  p  +  0  we  obtain 

<  BB*>  ■*  1,  hence  <  I  >  ->  I  . 

o 
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16.  HIGHER-ORDER  MOMENTS  <  IN> 


2 

Higher-order  statistics  start  with  <  I  >  .  At  the  time  of  writing,  controversy 
exists  about  the  behavior  of  <  >  as  a  function  of  the  parameters  k,  L,  and 

^  "(or  Cn^).  As  we  shall  see,  the  calculation  of  <  l2>  is  not  as  simple  as 
the  preceding  ones.  The  correlation  to  be  regarded  is  <  BB*  BB*> ,  and  it  is 
given  by  all  diagrams  with  4  axes  (2  full,  2  dashed)  that  we  can  form  by  insert¬ 
ing  beads  and  rungs  in  any  order.  Of  course  we  can  also  form  the  reduced 
diagrams  without  beads  by  summing  the  beads  to  give  a  factor  exp(-  4oL) ,  but 
first  we  do  not  do  this. 

2M 

Let  us  consider  a  diagram  with  M  features  (it  is  of  order  e  for  that  reason) , 
and  in  particular  let  us  concentrate  on  the  m-th  feature.  There  are  exactly 
10  possibilities  for  these  features  —  all  10  are  sketched  in  one  diagram  in 
Figure  3.  Our  notation  needs  some  further  development  in  order  to  write  down 
the  sum  of  ten  M-feature  diagrams  that  differ  only  in  the  m-th  feature  at  z^. 


1 

2 

3 

4 


Figure  3.  Ten  possible  rung  and  bead  features  at  any  location 
z  of  an  <  l2>  diagram. 


(i)  Notation  development:  We  will  designate  each  axis  by  an  index 
j  (j  -  1, 2,3,4).  The  B-axes  arc  1  and  3,  the  B*-axes  2  and  4  as 
indicated  in  Figure  3. 

(ii)  The  sum  of  all  Kp  vectors  from  p  =  m+1  to  p  =  M  that  connect  axis  j 

to  any  other  axis  will  be  denoted  by  QiUi*  Because  all  of  these  K 
vectors  match  each  other  two  by  two,  the  following  relationship  holds 
in  any  diagram: 


i(D 

V+i 


r(2)  *05 


;(4)  = 


0'“'  +  o'  -  O'  '  =  0 
m+1  ym+l  ha+1 


(16.1) 


(iii)  In  rder  to  write  down  the  sum  of  ten  factors  for  the  m-th  feature, 

T-Te  utilize  (13.5),  (13.7),  and  (13.8)  with  p  =  =0.  Common  to  all 

ten  factors  is 
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(16.2) 


2  2 

f  dz  /d2K  $(K )  x  (-  l)i+j+1  F  (ij) 
^g^2  ^7  vnj  mm  m 

o 


(lv) 

(V) 


where  Fm(ij)  stands  for  the  diverge  exponential  factors 
quoted  three  forms  for  the  m-th  factor. 

E^ch  of  the  ten  features  of  Figure  3  can  be  referred  to 
feature.  The  four  beads  occur  for  j  =  i,  the  six  rungs 
For  <  ij  >  =  <  12  >  ,  <  34  >  ,  <  14>  ,  <  23  >  we  obtain: 


in  the  above 

as  an  <  i j  > 
for  j  +  i. 


Fm(ij)  =  exp 


(L  -  zm)/k]  =  en(ij) 


(16.3) 


For  the  bead  contributions  j  =  i,  i.e.,  for  <  11>  ,  <  22>  ,  <  33>  ,  and  <  44  >  , 
we  obtain: 


F  (ii)  =  i 

m  4 1 


(16.4) 


For  <13>  and  <  24  >  we  obtain  the  somewhat  more  complicated  factors  Fm(13) 
and  Fm(24)  with 


F  (13)  =  e  (13)  exp  f-iK  2  (L  -  z  )/k] 
m  m  L  m  ml 

=  em*  (16.5) 

F  (24)  =  e  *(24)  e  * 
m  m  m 


where  we  have  further  simplified  the  notation.  With  these  conventions  we  can 
summarize  the  contribution  (16.?)  ior  the  m-th  factor  after  summing  over  all 
ten  diagrams  that  differ  in  the  m-th  factor  as  in  Figure  3,  to  obtain 


JcV 

16tt  2 


Cm  - 

13*1 


/  dz»j 


/ K  4>(K  )  x 
m  m 


(16.6) 


x  jfe  (12)  +  e  (34)  +  e  (14)  +  e  (23)-2l  -  \e  (13)e  +  e  *(24)  e  *l( 
Lm  '  m  m  m  J  [_  m  m  m  m  J( 


The  exponentials  m  (16.6)  can  be  converted  to  cosines,  which  are  abbreviated 
by  the  notations  (^(ij),  Cm  in  analogy  to  the  definitions  implicit  in  (16.3) 
and  (16.5).  Consider  (16.6)  in  conjunction  with  a  companion  set  of  10  diagrams 
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in  which  we  switch  the  roles  of  1  and  2,  3  and  4,  i.e.,  in  which  0^^  of  the 
new  set  equals  of  the  odd  set  and  vice-versa,  and  likewise  for  and 

-v(4)  ju 

We  ta^e  half  the  sum  of  these  two  sets  and  utilize  em(ij)  =  ^  (ji)  to 
order  the  indices.  The  result  is  that  the  { - }  factor  of  (16.6)  becomes 

{[c  (12)  +  C  (34)  +  e  (14)  +  e  (23)-2l  -  [c  (13)e  +  C  (24)e  *11 
(Lm  m  m  m  J  L  ®  m  m  m  J  f 

(16.7) 

Now  we  consider  a  third  set  of  diagrams  in  which  the  roles  of  1  and  4,  and  2 
and  3  are  switched.  Half  the  sum  of  this  set  and  (16.7)  yields  instead  of 
(16.6): 


tt 

16tt  2 


m-1 


/dz  f d^K  4>(K  )  x 

®  J  m  m 


j[c  (12)  +  C  (34)  +  C  (14)  +  C  (23) -2]  -  |~C  (13)  +  C  (24)]c  ( 

/Lm  m  m  m  JLm  m  jml 


(16.8) 


where  we  have  the  definitions. 


Cm  5  C0SK  <L  -  zm>/k] 

5  cos[v(^  -  (l  -  z„H 


(16.9) 


It  should  be  noted  that  we  need  the  other  M-1  factors  together  with  (16.8)  to 
describe  the  total  contribution  of  the^e  10  diagrams  differing  only  in  the 
m-th  feature.  The  reason  is  that  the  Qm+i  factors  can  contain  Kq  (q  >  m+1) 
vectors  which  make  (16.8)  also  dependent  upon  the  d^Kq  integration.  However, 
16.8)  is  still  not  correct.  To  see  this  let  us  consider  an  <  i j >  feature  at 
Zp  with  p  <  m.  It  depends  upon  Q^)  -  Q C j ) ,  which  is  written: 


Q(i)_q(j), 

P  P 


’  2  (k“>-  *‘J>)  +  ("m0-  km3’)  +  -  ^l)  (16'10) 


Unfortunately,  this  difference  varies^with  i,i  as  we  go  from  one  to  another 
Cm(ij)  in  (16.8)  becau^  it  contains  -  Km^).  Therefore  we  may  not  con¬ 

sider  the  p-th  factor,  for  p  <  m,  as  a  common  factor  preceding  all  terms  of 
(16.8);  it  changes  for  some  terms! 


Fortunately,  this  difficulty  can  be  circumvented  by  considering  the  sum  of  all 
(10)**  diagrams  obtained  by  allowing  each  and  every  one  of  the  M  features  to  be 
one  of  the  ten  in  Figure  3.  The  factor  for  m=l  will  be 


10 


-  (M-l) 


£  <[<^(12)  +  ^(34)  +  (^(14)  +  C1C23)-2]  - 

P1  ) 

-[^(13)  +  C1(24)]ci  i. 


(16.11) 


where  the  summation  indexed  by  P^  indicates  a  summation  over  all  10^"^  diagrams 
obtained  by  allowing  features  m  for  2  .s  m  <:  M  to  occupy  one  of  the  10  possibil¬ 
ities.  For  each  choice  of  the  second,  third,  ,  M-th  feature  there  is  a  sum 

{ - }  as  in  (16.11).  There  is  degeneracy,  i.e.,  not  all  of  the  10**-*-  factors 

{ - }  obtained  by  permuting  all  M-l  features  following  the  first  are  different. 

Now  consider  one  diagram  of  M  features  in  which  we  have  permuted  the  first 

feature.  We  obtain  one  of  the  terms  { - }  in  (16.11).  Now  permute  the  second 

feature.  As  we  have  shown  in  (16.10),  this  modifies  the  first  feature,  and 

consequently  we  obtain  (at  most,  barring  degeneracy)  six  different  terms  { - } 

of  (16.11)  to  go  with  each  second-feature  permutation.  However,  if  we  have 
the  entire  set  of  (16.11)  as  the  first  factor,  then  it  will  remain  unchanged. 

The  same  procedure  can  be  continued  for  m=2  up  to  m=M.  The  corollary  of  (16.11) 
for  the  m-th  feature  is 


10 


-  (M-m) 


£  [c  (12)  +  C  (34)  +  C  (14)  ->  C  ( 2 3 ) — 2]  - 
1 1.  m  m  m  m  J 


m 


Tc  (1J>  4  C  (24) |C  (  , 

L.  ic  m  j  m  s 


(16.12) 


where  Pm  symbolizes  one  out  of  the  10''  S’  diagrams  formed  by  permuting  features 

nri-1,  m+2, - ,  M,  Therefore,  the  contribution  tc  <l2>  for  the  sum  of  a] .1. 

M-th  order  diagrams  with  beads  and  rungs  is 


10 


-M(M-l)/2 


fi\  —fi\  n  ^4  f  dz  *«  ) 


m-l 


16c  ' 

o 


x  £  Jjc  (12)  +  C,  (34)  +  C  (14)  +  C  { 23 )  —  2 
ir*  j  j_  m  m  m  m 


p  p 
m 


(16.13) 


C  (13)  +  C  (24)1  C  , 
La  '  m  j  m  j 


who  re 
that 

not  t> 


one  must  count  all  permutations  P 
the  and  Cm  are  defined  In 

( j )  (  ^ 

ho-.-?  the  dependence  op.  0  ,  ,  -  0  4, 

K  n~*-l  m4  ]. 


m  including 
(16.9),  and 

explicit  Ly . 


the  degenerate  ones.  We  note 
that  the  notation  Cm(ij)  does 


There  is  nothing  very  special  in  (16.13)  about  the  restriction  to  <  I”>  dia¬ 
grams  aside  from  the  factor  10  and  the  number  of  Cm  factors  in  { - }.  When 

we  do  the  same  for  <  I^>  ve  have  the  following  choices  for  the  m-th  feature: 


(a)  2N  beads 

fl>)  N(N-1)  rungs  of  (13.8)  type, 

(c)  rungs  of  (13.7)  type. 

The  contribution  to  <  I^>  by  all  M-th  order  diagrams  is  given  by  (16.13) 
provided  we  replace 

10-«(M-l)/2  by  [N(2H+  1)rM<M-l)/2_ 

C  (12)  +  —  +  C  (23) j  by  N  terms  (one  for  each  BB  rung) 
m  m  j 

-2  by  -  N  (a  term  -  1/2  for  each  bead),  \  (16.14) 

[c  (13)  +  C  (24)1  by  N(N-l)  terms  (all  BB  and  all  B*B*  rungs). 

L  m  m  j 


Although  (16.13)  and  (16.14)  are  rather  complex,  it  is  possible  to  draw  some 

preliminary  conclusions  from  them.  To  do  so,  we  first  examine  the  { - }  factor 

because  it  apparently  acts  as  a  filter  function  upon  <J> (K) .  We  shall  show  that 
it  is  indeed  a  "high-pass"  filter. 

Consider  one  { - }  factor  of  (16.13),  to  be  specific,  and  expand  the  cosines 

into  power  series  of  K^.  it  can  oe  seen  by  utilizing  (16.1)  that  the  terms 
proportional  to  Kffi2  cancel  each  other  (as  do  the  terms).  In  fact,  we  can 
prove  this  for  (16.14).  We  note  immediately  that  the  unity  terms  of  the  cosines 
cancel.  The  terms  come  from  the  second  terms  of  Cm(ij)  setting  Cm  =  1. 

They  can  be  written  (omitting  constant  factors)  as: 


where  the  last  form  is  the  generalization  oi.  flh,])  ana  therefoie  zero.  As  a 
consequence,  reduces  to 


5! 


(16.15) 


{ - }  -  K4  (L  -z  )2/k  +  K4  AQ4..  (L  -  2  )4/k4  + - , 

mm  m  nrhL  m 

where  the  next  term  is  of  0 (Km) ^ ,  etc.  We  have  abbreviated  the  sum  over  all 

i  $  j  of  the  4-th  powers  of  -  QmH  by  AQ^.j.  Now  note  that  the  2nd  term 

-*■4  2  2 

contains  an  extra  factor  AQm+i  (L  -Zjq)  /k  above  the  first.  If  this  factor  is 
much  less  than  unity,  then  -  clearly  -  the  first  term  of  (16.15)  is  the  leading 
term  of  { - }.  Let  us  choose  a  A  «  1.  It  follows  that 

A(k/L)1/2 

dKK  4»(K)  { - }  « 

o  o 


oo 

/ 


dKK  *(K)  { - }, 


(16.16) 


consequently  we  note  that  { - }  filters  out  all  wavenumbers  ouch  less  than 

(k/L)*/2.  This  is  rather  useful  in  optics  where  (k/L)l/2  and  as  a 

consequence  we  note  that  { - }  in  optics  is  estimated  roughly,  and  parametric¬ 

ally  by 


dKK  4>(K) 


*  k7/6  l11/6c  2 

n 


(16.17) 


where  we  have  utilized  <t(K)  «  L0“2^3  This,  too,  is  a  well-known  result 

noted  by  diverse  Russian  authors[l4].  In  radiowave  propagation,  (k/L)l/2  < 
L0~l,  hence  the  matter  is  academic  there  because  the  left-hand  side  of  (16.17) 
is  then  hardly  different  from  a  L. 

4  2  2 

The  second  term  of  (16.15)  is  smaller  than  the  first  when  ^AQ^^  L/k  >is 
much  less  than  unity.  Equivalently,  we  can  set  kL  <  Qm2/k2>  much  less  than 
unity.  We  utilize  <Qm2/k2>  ^  a2<(60)  '>,  as  in  Section  8,  to  find 


<Q2  L/k  >  *  kL2  k  1/3  L  ~2/3  e  2 
in  mo 


(16.18) 


Note  that  this  is  an  overestimate.  Jn  applying  (14.7),  we  should  now  -  as  in 
(16.17)  -  have  replaced  the  lower  bound  of  the  "averaging"  integral  by  a  wave- 
number  Kq  'v  (k/L)i/2t  Nevertheless,  if  we  do  so,  we  still  obtain  (16.18)  be¬ 
cause  the  Integral  of  K^  $(K)  is  determined  by  the  immediate  vicinity  of  K  =  Km, 
and  therefore  only  weakly  dependent  on  the  lower-wavenumber  cutoff  'v  (k/L)l/2. 
Note  that  by  regrouping  the  terms,  replacing  L0~2/3  e2  by  an(j  utilizing 
at  2  ^  k7/6  Lll/S  Cn2,  we  can  rewrite  (16,18): 


^  ^2  t  /i  ^  „  /  2  .  .  1/  6  2 

<  Q  L/k  >  s,  (k  L,  k)  oe 
m  m 


(16.19) 
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In  optics  (<m  L/k)^  is  of  order  unity.  For  example,  for  k  *  10  m  , 

L  *  10  km,  and  iQ  •  6  mm,  one  obtains  (icm2  L/k)^/^  ^  3,  clearly  tbe  size  of 
the  second  term  in  (16.15),  hence  the  validity  of  (16.17),  is  strongly  deter¬ 
mined  by  the  parameter  oe2.  I.e. ,  (16.17)  is  self-consistently  small. 


17.  THE  HODIFIED-RYTOV  APPROXIMATION 


The  name  of  "Rytov  approximation"  appears  to  be  tied  rather  loosely  to  a 
number  of  concepts.  It  is  associated,  for  instance,  with  the  first  approxi¬ 
mation  of  the  Russian  "method  of  smooth  perturbations  (MSP)"[10,31] .  In  this 

case  one  takes  Equation  (2.3),  substitutes  E=E  exp  ip  to  obtain, 

0 


W  +  (W2  +  2ik*V*  +  k2  6e 


0, 


and  then  the 

h 


2 

(ViJ>)  term  is  dropped.  The  result  is 
=  f  dz  fih  e  -«2a-2)/2k 


j'K*p 


(17.1) 


(17.2) 


The  statistics  of  B  *  exp  i are  easily  computed.  In  particular  let 
!|»1  *  Xi  +  <(>x*  It  then  follows  that 


<  IN> 


.N 


2Nxi 

<  e  >  = 


1 1  exp[2N 2<xl> 


2  2  7/6  11/6 

a  =  0.307  C  k  L 
e  n 


(17.3) 


These  results  are  quite  well  known,  and  the  derivations  are  relatively  simple 
and  discussed  elsewhere [10, 11, 15, 16] .  It  is  quite  obvious  that  (17.2)  can¬ 
not  be  correct,  as  is,  because  energy  is  not  conserved:  <I>  =  ID  exp(2a^). 
Also,  it  is  well  known  at  this  date  that  the  normalized  irradiance  variance 
o£  defined  by 

=  [  <I2>-  <1>  2  ]/  <  I>  2,  (17.4) 

2 

given  in  (17.3)  explicitly  as  exp(4o£)  -  1,  or  that  the  log-amplitude  variance 
<xf^>=  does  not  behave  as  measured  when  approaches  and  exceeds  unity[15]. 

However,  the  formalism  of  Section  16  contains  a  version  of  the  Rytov  approxi¬ 
mation  that  does  conserve  energy.  Let  us  reconsider  (16.13)  and  (16.14)  for 
the  case  that  we  may  replace  in  (13.1)  by 


Fn(VWZm>  =  exp[-iKm(L-zm)/2k] 


(17.5) 


That  is  to  say, 
equal  to  unity. 
Cm(ij)=l.  This 
bution  of  every 
(16.13)  reduces 


we  ignore  the  exponential  exp  [-iKj,,  •  Qm-i (1-%) /k]  by  setting  it 
If  we  then  substitute  (17.5)  into  (16.2)  -  (16.8)  we  obtain 
constitutes  an  enormous  simplification  cf  (16.12).  The  contri- 
term  in  the  sum  is  identical,  and  equal  to  2(1-Cm).  Hence, 
to 
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(17.6) 


M 

n 

m=l 


m-1 


2  2 
k  e 

4w 


/  «./ 


dK  K  4>(K  ) 
mam 


cos  [K 


m 


•(L-z  )/k] 
m 


}■ 


2  H 

(4op  /M! 


where  we  have  utilized  (7.1),  and  the  well-known  relationship  between  the 
above  factors  and  the  calculation  oc  a2  in  (17.3).  Because  v°  have  demonstrated 
that  <  I2>is  given  by  summing  (17.6)  from  M=0  through  all  positive  integer 
values,  we  note  that  we  obtain 


<  I '>  =  I  exp[4oc  '] 


(17.7) 


*N, 


Similarly,  for  <  I  >  we  note  by  looking  at  (16.14),  that  { - }  is  replaced  by 

N(N  -  1)  (1  -  Cm)  •  In  this  case  the  M-th  term  of  <I^>  is  given  by  (17.7) 
with  one  difference:  each  factor  k2^/^  *oe  replaced  by  N(N  -  l)k2e2/3Tr. 

This  yields  [2N(N  -  1)  oe2]  ^/H!  as  the  contribution  of  til  M-feature  diagrams 
to  <  T^>  .  Therefore, 


<  IN>  =  I/  exp  [2N(N  -  1) 


(17.8) 


Note  that  energy  is  conserved.  Note  also  how  closely  (17.8)  resembles  (17.3). 
In  fact  if  <Pl,  in  (17.2)  is  replaced  by  4<1  -  oe2>  the  statistics  of  (17.3)  be¬ 
come  identical  to  (17.8).  Hence,  the  approximation  (17.5)  yields  the  modified- 
Rytov  approximation  in  optics: 


E 


E 

o 


exp  (^  -  oe  ), 


(17.9) 


where  4^  is  given  in  (17.2). 

Under  what  conditions  is  (17.8)  or  (17.9)  valid?  No  really  convincing  answer 
appears  to  have  been  given  to  date,  other  than  that  oe 2  must  be  small.  However, 
it  appears  to  us  that  a  good  criterion  can  be^founcj  by  examining  when  (17.5)  is 
a  good  approximation.  I.e.,  when  does  exp[-iKm  -  (L  -  zm)/k]  behave  as 

unity  in  the  formalism?  The  answer  is  very  easily  estimated  by  retracing  the 
development  of  (16.8).  We  simply  examine  [instead  of  (16.8)], 


,2  2 

ke 


16tt 


L 

J '  dz  jd^K  $00  jl  ~  cos  [K^  (L  -  z) /k]j  exp [-iK* AQL/k] 
o 


(17.10) 


because  this  will  overestimate  the  influence  of  the  last  exponential  factor  if 
it  affects  the  integral.  We  convert  (17.10)  into 
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(17.11) 


As  noted  in  developing  (16.17),  the  influence  of  the  filter  factor  { - }  is 

simply  to  eliminate  contributions  to  (17.11)  from  K  <<  (k/L)l/2,  Hence  (17.11) 
is  estimated  by  setting  K  $(K)  =  15.7  (because  the  lower  bound 

'v*  (k/L)l/2  far  exceeds  che  wavenumber  where  KL0  =  1)  to  obtain 


15.7  C  \k2 
_ n 

Sir 


f  dK  K_8/3 

/t  \l/2 


J  (K  AQL/k) . 
o 


^  c  2  k7/6  l11/6  /-  x"11/6  (k/L)1/2  A  QL/k] 

1 


(17.12) 


2 

We  obtain  oe  if  the  Bessel  function  does  not  contribute  appreciably  to  the 
integral.  Therefore  a  sufficient  condition  is  (k/L)l/2  AQL/k  <<  1,  or  „ 

(AQ/k)2  kL  <<  1.  However,  we  have  already  estimated  < (AQ/k) 2 >  'v  2L  Cn 

in  (14.7).  Consequently,  the  criterion  for  modified-Rytov  (17.9)  is 

L2  <m1'\~213  k  e2  «  1,  or  (k2  L/k)1/6  a  2  «  1  (17.13) 


which  is  identical  to  (16.19).  Note  therefore  that  this  development  confirms 
(16.17)  which  we  have  now  rederived  more  carefully.  Note  also,  as  explained  in 
connection  with  (16.19),  that  (icm2  L/k)l/6  ^  3  or  less  so  that  the  condition 
(17.13)  is  not  very  different  from  the  well-known  condition  o£2  << 

1.  We  note 

that  (17.13)  relaxes  this  condition  slightly.  Finally,  (17.13)  is  also  equal 
to  the  condition  (8.10)  for  the  Moli&re  approximation,  although  this  latter 
approximation  is  more  severe  in  the  multiple  -  scattering  regime  (see  Section 
19).  Our  interpretation  is  that  (Km2  L/k)l/6  ae2  can  get  somewhat  closer  to 
unity  for  the  modified-Rytov  approximation. 
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18.  THE  C  -  0  AND  C  (ii)  »  1  APPROXIMATION 
m  m 


When  Cm  =  0  and  Cm(ij)  *»  1,  and  all  diagrams  for  which  that  is  not  the  case 

are  ignored,  we  obtain  a  result  that  appeared  previously  to  us  to  be  limited 

to  low-frequency  and  long-distance  propagation,  hence  to  the  radiowave  case. 

We  intenJ  to  examine  it  more  carefully.  First  of  all,  we  define  what  we  mean 
by  "C*  -  0". 

2 

When  L  »  kLo  ,  it  follows  from  the  comments  following  (16.17)  that  the  factor 
Cm  plays  no  role  of  importance.  To  state  this  another  way,  when  all  of  the 

Cjn(ij)  *  1,  then  the  filter  factor  { - }  will  weight  K  4>(K)  as 


f  dKK  4>(K)  [1  -  cos  (K^/k)], 
o 


(18.1) 


if  we  replace  L  -  z  by  L  to  make  matters  more  extreme.  Tatarski  [10,  Section 
7.4)  has  discussed  the  effect  of  filter  factors  such  as  1  -  cos  K2L/k.  It  is 
easily  seen  that  the  second  term  of  (18.1)  is  negligible  when  L  >>  kL02.  Just 
consider  $  as  a  function  of  k2l02  **  x,  so  that  (18.1)  can  be  written  (aside 
from  constants) , 


/ 


dx  * t  (x) 


[1  -  cos  (x  L/kl^2) ] 


(18.2) 


where  we  have  used  the  notation  $(10  =  $'(x)  [which  is  proportional  to 
(1  +  x)~H/6],  The  first  term  of  (18.2)  is  essentially  the  one-dimensional 
covariance  of  6e,  C(0),  at  the  zero  of  its  argument.  The  second  term  is  essen¬ 
tially  C  (L/kLQ2).  Consequently,  when  all  Cm(ij)  =  1,  we  may  set  the  filter 
factor  (1  -  Cm)  »  1,  and  the  error  is  of  order  C  (L/kLQ2)  /C  (0).  In  that 
sense,  we  say  "Cm  ■  0."  Note  that  the  variable  that  serves  as  argument  of  C 
has  the  dimension  of  y/L0,  where  y  is  *  separation  distance  corresponding  to 
the  correlation  <  <$e  (r  +  y)  <5e  (r)>  averaged  over  all  r.  When  L  >>  kLo  ,  the 
above  rat*u  is  extremely  small.  Therefore  <  13 >  and  <  24 >  rungs  (which  are 
responsible  for  C^)  play  hardly  any  role  in  this  case. 

However,  instead  A  working  with  (16.13),  and  (16. ]4)  we  shall  use  a  slightly 
different  formulation.  We  will  consider  <  |BB*|  >  as  the  sum  of  all  possible 

diagrams  reduced  by  summing  out  beads  as  described  in  Section  14.  Instead  of 
forming  (16.13)  we  form 


-4aL 
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m=l 
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16tt 
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m-1 


dz  4>(K  ) 
m  m 


F  (ii) 
m 

(18.3) 
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as  the  sum  of  all  M-rung  diagram  contributions  to  <BB*BB*>,  where  Fm(ij)  = 
Cm(ij)or  Fm(ij)  =  Cm(ij)  C^,  analogous  to  (16.3)  and  (16.5).  Thus  for  one 

diagram  we  obtain  a  product  Fi(ii,jj)  F2(i2»j2) - ln  tue  inte“ 

grand  of  (18.3)  where  im,jm  refers  to  the  K  vectors  of  Qm_i  between  axes  im,  jm. 
I.e.,  1  <  im,  jm  <  4,  and  jm  ^  im,  because  beads  have  been  summed  out.  The  sum¬ 
mation  sign  with  index  Pjj  indicates  the  sum  over  all  topologically  different 
diagrams.  This  is  in  contrast  to  (16.13)  where  we  summed  some  diagrams  more 
than  once!  Here,  there  is  no  degeneracy  factor. 

The  formulation  (18.3)  serves  to  discuss  the  elimination  of  contributions,  for 
which  Cm(ij)  £  1.  In  the  above  discussion  of  (18.2)  we  showed  that1  any  term  of 
(18.3)  containing  a  factor  Cm(13)  and/or  Cm(24)  in  the  integrand  is  small. 

I.e.,  it  is  of  order  [C(L/kLo2)/C(0) ]n  if  there  are  n  sucl.  such  factors,  com¬ 
pared  to  what  we  obtain  if  all  C^Uj)  =  1. 


Now  we  estimate  the  contributions  of  Cm  (14),  Cm  (23),  C^,  (.12),  C-,  (34)  when 

^(j) 

4-  Qm  for  these  four  factors.  An  overestimate  is  given  by  considering 


*(K)  exp  [iK  -AO., 
m  m  nH-1 


(L  -  z  )/k] 
m 


(18.4) 


as  an  entirely  independent  factor.  Other  Km-dependent  factors  are  also  oscil¬ 
latory,  and^their  effect  is  estimated  anyway  by  replacing  L  -  zm  by  L  and  by 
modifying  A  Qm+i  to  include  the  ether  Independent  terms  in  the  exponentials. 
Hence  we  compute 


*(K)  exp  [iK-AQ  L/k]  =  2tt  /  dKK  $(K>  J  (KAQ  L/k) 

I  Q 


”/ 


(18.3) 


for  AQ  ^  0  and  compare  it  to  what  we  obtain  for  AQ  -  0  to  estimate  the  relative 
contribution  of  compared  to  unity.  For  the  Kolmogorov  spectrum  (10.10) 

we  can  rewrite  (18.)),  except  for  inessential  constants, 


00 

/dx  (1  +  x)"11/6  Jq  (yxfx)  =  4  yL5/t  K5/6  (yx) /211/6r  (11/6)  (18.6) 

o 

where  x  =  K^LQ2  and  y^  =  AQ  L/kL0.  Thus,  the  ratio  of  Cm  (AQ)  /Cm  (0)  is  esti¬ 
mated  by  the  ratio  of  (18.6)  for  y1  0  to  (18.6)  for  y^  =  0,  i.e., 


C  (AO)  . 
m  _____  6__ 

C  (0)  %  11 

m 


i/6  5/6 


K 


5/6 


(7)) 


(18.7) 
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Here,  as  in  (18.6),  we  use  the  NBS[13]  notation  for  Bessel  and  ganma  functions. 
Furthermore  (^(AQ)  is  a  convenient  notation  for  the  Cm(ij)  under  consideration 
here.  Now,  referring  to  (14.7)  we  estimate  p  by  estimating  the  rms  of  AQ,  i.e., 
we  know  that  the  variance  of  in  the  sense  discussed  in  Section  6  is. 
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9  9  9  9 
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■  (K  L  )1/3  (aL)  (L/k  L 
mo  o 


K  1/3 
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l-8/3.j 
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(18.8) 


2 

Clearly  <  y^>  >>  1  in  the  radiowave  regime  under  consideration,  unless  aL  <<  1 
(in  which  case  we  have  well-understood  single  scattering) .  Therefore,  except 
for  a  numerical  factor,  (18.7)  reduces  to 


C  (AQ) 
m _ 

C  <°> 
m 


ry-i 


(18.9) 


Let  us  return  to  (18.2)  and  actually  do  the  estimate  for  the  Kolmogorov  spec¬ 
trum.  We  obtain 


on 

C  ^  f dx  (1  +  x)  lly^6  cos  xv0 
m  J  £■ 

O 

y2  =  L/kLo2, 


(18.10) 


for  the  effect  of  the  cosine  term  in  (18.2)  divided  by  the.  first  term.  Using 
the  complete  gamma-function  notation  for  imaginary  exponentials,  see  Reference 
[13]  (6.5.7  and  6.5.8),  we  convert  (18.10)  into 

C  'v-  y,5^6  (cos  y0  C  (y  ,  -  5/6)  +  sin  y  S  (y  ,  -  5/6)] 

Z  1  1  1  1  (18.11) 

*  (1/kL  2)“2 
o 


"vnere  we  have  used  asymptotic  forms  of  the  C  and  S  functions  for  large  y?  = 

L/kL02. 

We  shall  postpone  discussion  of  the  cumulative  effects  of  errors  (18.9)  and 
(18.11).  First,  we  treat  two  remaining  problems:  (a)  the  contribution  of  any 
M-rung  reduced  diagram  for  which  =  0  and  Cm  (ij)  =  1,  and  (b)  the  bookkeeping 
of  all  possible  diagrams.  The  first  problem  is  practically  trivial  when  deduceu 
from  (18.3),  because  Fm  (ij)  =  i  under  the  assumptions  just  stated,  and 
(-  l)i+j+l  =  1.  Hence  we  obtain 


e 


(18.12) 
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f 

dKK  4>(K)  = 
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4CL  (2aIi)M 
M! 


as  the  contribution  of  an  M-rung  reduced  diagram.  The  bookkeeping  problem  is 
as  follows:  First,  we  note  for  <I^>  that  Cm  =  0  means  no  <  13>  and  no<24> 
rungs  are  permitted.  Then,  Cm(ij)  =  1  implies  that  there  are  only  one  of  the 
two  following  possibilities  in  any  single  diagram. 

(i)  Only  <12>  and/or  <  34>  rungs. 

(ii)  Only  <  14 >  and/or  <  23 >  rungs. 

Any  mixing  will  cause  at  least  one  Cm(ij)  ^  1.  In  order  to  avoid  counting  a 
diagram  more  than  once,  we  add  up  contributions  (18.12)  as  follows: 

The  no-rung  diagram:  The  contribution  of  all  beads  summed  out  to  yield  a 
no-rung  reduced  diagram  is  obviously 


exp  (-  4aL)  (18.13) 

One  pair  of  axes  without  rungs:  In  this  case  we  assume  rungs  only  of  <  l4>or 
of  <  23 >  or  of  <  12 >  or  of  <  34 >  type.  In  each  of  these  four  cases  we  add 
up  (18.12)  for  M  =  1  to  M  •  °°,  to  obtain  for  all  four. 


<4  exp  (-  4aL)  [exp  (2aL)-l]  (18.14) 

A  double  set  of  rungs:  In  this  case  we  assume  that  there  are  rungs  either  on 

<  12 >  and  <  34 >  at  the  same  time,  or  on  <  14>  and  <  23>  .  In  other  words, 
both  of  the  above  cases  are  excluded.  Because  all  possible  relative  orders 
between,  say,  <  12>  and  <34>  rungs  occur,  we  may  separate  (18.12)  into  a 
double  sum  over  mi2  and  m34  where  M  =  mi2  +  m34  (mij  Is  the  number  of  rungs 
between  axes  i  and  j) .  We  permute  all  M  rungs  but  exclude  permutations  of 

<  12 >  and  of  <  34 >  rungs  among  themselves  in  order  to  generate  all  M-rung 
diagrams  of  this  double-rung  category.  Hence  we  obtain: 


2  exp  (-4aL)  V 
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M 


=  2exp  (-4aL)  [exp  (2aL)-l]  . 


(18.15) 


The  factor  2  is  present  in  (18.15)  so  as  to  include  the  <14>  +  <23>  group  .  We 
note  that  <  l2>  / < 1 >  2  is  therefore  just  the  sum  of  (18.13),  (18.14),  and  (18.15) 
since  we  have  exhausted  all  possible  diagrams  and  carefully  avoided  counting  any 
contribution  more  than  once.  The  result: 

<  1^>  =  <1>A  [2  -  e^p  (-4aL)]  (18.16) 
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N 

Now  we  do  the  same  procedure  for  <  1  >  .  Up  and  through  (xS.12)  there  is  no 

change  except  exp  (-  4aL)  becomes  exp  (~  2NaL) .  However,  the  bookkeeping  prob¬ 
lem  changes.  Because  =  0,  we  have  rungs  only  between  B  and  B*  axes.  Let 

us  try  to  add  up  (18.12)  for  all  M-rung  diagrams  where  n  pair  of  B  and  B*  axes 

have  no  rungs  [obviously  we  mean  reduced  diagrams  so  that  the)  is  a  common 
factor  exp  (2NaL)].  We  may  not  have  any  mixing,  i.e.,  we  cannot  have  two  rungs 
connecting  one  B  to  tuo  B*  axes,  or  vice  versa,  otherwise  at  least  one  C,n(ij)  ^  1. 
For  one  choice  of  n  pairs  of  conjugate  B  and  B*  axes  there  are  (N-n) !  ways  of 
pairing  the  remaining  N-n  pairs  of  B  and  B*  axes  nwo  by  two.  Let  us  assume  for 

one  choice  of  pairing  that  there  are  mi,  m2,  -  nn$_n  rungs  on  the  N-n  pairs  we 

have  chosen.  Obviously  mi  +  -  +  mn-.n  =  M.  Then,  analogous  to  (18.15)  we 

obtain 


(N-n) ! 


-2NaL 
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00 


-MaL 


(N-n) !  e2NaL  (l  - 


-2ctL\N-n 
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(18.17) 


for  all  M-rung  diagrams  that  leave  a  single  choice  of  n  pairs  of  3,B*  axes  un¬ 
connected  to  each  other  by  rungs.  We  need  one  more  numerical  factor  in  (18.17) 
before  summing  over  n  from  n  *  0  to  n  =  S,  namely  the  number  of  ways  that  we 
can  choose  n  pairs  B  and  B*  from  N  pairs  B  and  B*.  This  number  is  easily  seen 
to  be 


N2  (N-l) 2 (N  -  n  +  l)2 

~2  ~2  2 

1  2  - n 


_ (Nil! _ 

(n!)2  [(N-n)!]2 


(18-10) 


N  N 

Thus,  we  obtain  <  I  >  /  <I>  by  summing  the  product  of  (18.17)  with  (3  8.18) 
over  n  from  n  =  0  to  n  =  N: 


j.N  >  =<I>N 
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(18.19) 


Note  that  (18.16)  is  indeed  the  result  (18.19)  for  N  =  2.  Note  also  that  N  =  1 
yields  an  identity.  We  have  already  derived  <  I>  =  I0,  hence  <  I>N  can  be 
replaced  by  IQN  in  (18.19). 


This  result  is  extremely  interesting  because  it  can  be  shown  that  (18.19)  _or- 
responds  to  the  Rice  distribution  for  I  defined  by  the  parameters  10,  and 
1^  =  I0  exp  (-  2aL) ,  jith  orobability  density 


o(I)  =  (I  -  1  ) 

O  L 


exp 


l~(T  + 
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(18.20) 
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First  use  NBS  formula  22.3.9[13]  to  show  that  (18.19)  is  essentially  a  Laguerre 
polynomial,  then  we  compute  the  integral  of  IN  p  (I)  for  I  =  0  to  I  =  00  in 

(18.20)  by  using  NBS  formula  22.10.14  to  obtain  agreement  with  (18.19). 

Another  equivalent  way  of  stating  either  (18.19)  or  (18.20)  is  to  say  that  the 
electric  field  E  is  the  sum  of  a  constant  part  El  =  Eo  exp  (-aL)  and  a  Rayleigh- 
distributed  part  6E  with  <6E>  =  0  and  <|6S|2>  =  i0  -  II* 

N 

Thus,  we  note  that  the  selective  summation  of  diagrams  contributing  to  <1  >  , 
ignoring  those  for  which  Cn  (ij)  t  1  occurs  (as  well  as  by  setting  Cm  =  0) , 
yields  the  Rice  distribution!  Now  we  wish  to  establish  conditions  under  which 

(18.20)  holds.  We  utilize  (18.9)  and  (18.11)  but  now  we  estimate  the  cumulative 
effects  of  the  errors.  First  of  all,  we  note  that  the  role  of  o2  (as  defined 
in  Section  4)  is  played,  again,  by  aL-  So  we  will  set  cj2  =  aL  and  follow  the 
averaging  procedure  of  Section  8  where  needed. 

Consider  first,  the  effect  of  <BB>  and  <B*B*>  rungs,  e.g-,  of  <13>  and 
<24>  in  <l2>  .  Each  of  these  gives  rise  to  a  factor  -  Cjn  (ij)  Cm,  as  noted 
in  (18.3).  The  error  m  each  is  given  by  (18.9)  unless  Cm  (ij)  =  1.  However, 
the  only  way  that  Cm  (ij)  can  be  unity  is  that  the  rung  in  question  is  tne  last 
one  in  which  case  it  generates  as  many  factors  Cp  (i'j')  /  1 (Cp  <  m)  as  tnere 
are  rungs  ending  at  zp  on  either  i  or  j  (i1  and  j'  are  other  index  pairs,  not 
necessarily  the  same  as  i,j).  The  only  way  in  which  an  error  effect  (18.11) 

can  be  isolated  from  08.9)  is  when  any  pair  of  two  B  or  two  B*  axes  are  con¬ 
nected  by  one  and  only  one  rung,  and  neither  of  the  axes  of  the  pair  is  con¬ 
nected  to  any  third  rung.  For  <I^>  ,  the  error  appears  to  be 


<IN_2> 


,  2  2 
k  e 

8ir 


j  dz  J  dKK  *(K)  C1 
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N-2 

compared  to  <1  ^  if  we  allow  one  isolated  <BB>  rung.  The  error  can  also 

be  (L/kL02)-4  <jN-2>  if  we  have  .a  <BB>  and  a  <B*B*>  rung  isolated  from 
the  2N  -  4  others;  however  this  error  is  negligible.,  Therefore,  the  cumulative 
error  appears  to  be  0  (L/kL0^)~2. 

Now  we  consider  the  effect  of  Cm  (1?\  Cm  (34),  (14),  and/or  Cm  (23)  not 

equal  to  unity  [or,  for  <I^>  ,  the  effect  of  Cm  (ij)  /  0  for  any  of  the  ij  in 
the  second  row  of  (16.14)].  It  is  best  examined  by  returning  to  (16.13)  and 
(16.14)  and  reexamining  the  filter  factor.  Let  us  take  a  diagram  for  which  all 

Qii+1  "  QmJl  =  0.  By  placing  the  m-th  rung  at  all  possible  BB*  connections,  we 
generate  the  N2  Cm  (ij)  of  the  seconJ  li«c  of  (6.14),  and  of  course  all 
Cm  (ij)  =  1-  However,  for  each  B  v,N  at  most),  we  can  generate  N  -  1  connections 
to  E*  axes  ’hat  will  yield  a  Cm_f  (ij)  f  1,  and  only  one  that  yields  Cm_]_  (ij)  =  1. 
Hence,  in  the  (,m-l) st  filter  factor  there  are  at  worst  a  ratio  of  1/(N-1)  of 
factors  (ij)  /  1.  If  the  features  m+1,  m+2,  - ,  M  are  not  so  chat 
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Qm+1  “  Qmil  =  0,  we  cannot  get  as  large  an  error.  Hence  the  error  in  <IN> 
appears  to  be 


0  (Ny1/3  e'y)  and  0  ^L/kLQ2) 

.  2.  n  T  3  1/3  T  -8/3  e  2 

<y  >  'v  L  k  L  e 

J  mo 


(18.22) 


Therefore,  (18.20)  is  a  good  approximation  for  all  applications  in  which  only 
lower  order  moments  <1^>  with  N  <<  exp(y)  are  involved.  When  y  »  1,  this 
is  not  a  serious  restriction.  Note  that  L  >  kLQ2  implies  <y2>  >>  1,  hence 
(18.21)  suffices  for  (18.20). 


63 


19.  PARAMETER  REGIMES:  A  GRAPHIC  REPRESENTATION 


So  far,  we  have  derived  a  number  of  results,  each  under  certain  assumptions 
and  therefore  restricted  by  errors  of  the  type  given  in  (18.22).  Such  errors 
define  relationships  between  end  L  (namely  by  setting  the  error  parameter 
equal  to  unity) .  We  can  make  matters  much  clearer  b>  introducing  a  graphic 
plot  of  vs.  L  as  in  Figure  4.  The  actual  plot  has  e*  as  ordinate,  and  L/L< 
as  abcissa.  We  plot  each  coordinate  logarithm: cally  in  units  log  kL0  so  that 


L/L0 


Figure  4.  Turbulence  strength  vs.  pathleng:h. 


power  laws  are  straight  lines.  Note  that  the  plot  contains  regions  labeled 
I,  II,  III,  IV  (some  with  indices)  's  well  as  regions  labeled  otherwise.  All 
regions  are  bounded  by  lines.  We  discuss  the  general  regions  first: 

The  diffusion  region  Aq  £  I .  In  the  upper  righthand  corner  of  Figure  4,  we  have 
plotted  the  line  Aq  =  1,  using  (14. 7). _  Whenever  Aq  £  1,  we  no  longer  have 
n  <(60)2>  much  less  than  unity,  nor  Qm^/k2  <<  1.  Hence,  diffusion  occurs  be¬ 
cause  n  is  the  number  of  scatterings  and  6$  =  H - f-  6^n  is  an  inconerent 

sum  with  zero  mean  and  variance  n  <(60)2>.  The  transition  from  (5.4)  to  (6.10) 
requires  Aq  «  1.  Therefore,  all  results  in  this  report  require  exclusion  of 
this  region. 


The  multiple -olu8 ter  region:  A^e  £  1.  This  is  the  region  lying  above  the 
horizontal  line  e2  =  (kL0)-  .  It  corresponds  to  A be  K  15 it  is  no  longer  per¬ 
missible  to  regard  the  n-point  correlation  of  6e  as  a  product  of  binary  cor¬ 
relations.  All  of  our  statistical  results  exclude  Agc  £  1. 


The  sagittal-error  region  TVgg:  As  K  1.  For  the  single-scattering  regime  it  is 
defined  by  L  £  k^2.0^  in  optics  (not  denoted  in  Figure  4  because  it  does  not  lie 
in  the  optics  region)  and  L  £  for  radiowaves  (a  vertical  line) .  For  the 

multiple-scattering  regime  As  is  given  by  (14.8).  It  is  sketched  as  the  border 
between  IVg  and  IV«;g.  In  region  IVgg.  the  exponential  function  in  (6.10)  is  no 
longer  a  good  approximation,  and  consequently  most  (but  not  all)  of  the  results 
we  derive  are  not  vaiid  in  IVqg.  The  sagittal  error  does  not  affect  optics 
results.  Nov  we  shall  discuss  the  significance  of  the  major  dividing  lines. 


The  optics-radiowaves  boundary:  L  =  kL0  .  We  have  noted  that  in  optics  we  usually 


t 

have  k£02  <  L  <<  kLo  because  k 


10*  m"l 


1  is  a  very  high  wavenumber  (kL0 

In  radiowave  propagation  we  have  k  £_10  m  ,  but  L0  £  10  m.  Although  it  is  pos- 


2  >  10**  km)  . 


-  *  W  v.  '  V  u  1 

sible  -  easily  so  -  to  have  L  <<  kL0z  in  this  case,  it  is  also  true  that  a  L  «  1 
in  that  case  (see  later).  The  really  interesting  results  occur  for  L  >>  kLr2 
which  -  apparently  -  can  only  occur  at  lower  frequencies,  hence  in  the  radiowave 
regime.  The  L  =  kL02  line  thus  defines  two  optical  regions  I  and  III,  and  two 
radiowave  regions  II  and  IV. 


Tas  single-multiple  scattering  boundary  a L  -  1.  When  aL  «  1,  we  note  that 
(18.20)  reduces  to  a  delta  function  6(I-I0),  and<B>»l  in  (14.6).  It  is  well 
known  that  the  first-order  corrections  (although  they  have  no  physical  meaning 
because  energy  is  not  conserved)  are  due  to  the  approximation  B  =  1  +  in 
(3.1).  This  therefore  is  sufficient  for  defining  single  scatter-'ng.  Consequently, 
aL  %  1  must  include  multiple  scattering  effects.  Hence  aL  =  1  [a  is  defined  in 
(18.12)]  delineates  the  regions  of  multiple  scattering  (III  and  IV),  and  the 
regions  I  and  II  are  therefore  single-scattering  regimes. 


It  is  now  much  easier  to  skim  through  the  results  derived  so  far  and  to  delineate 
their  regions  of  validity: 

(i)  The  WKB  approximation  derived  in  Section  (7)  is  limited  by  (8.10) 

in  the  multiple-scattering  regime,  and  by  (7.5)  in  the  single-scatter¬ 
ing  regime.  Hence  it  is  valid  in  I  and  in  IIIW  (the  left  part). 
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(ii)  T'ivi  modi  fie  d-Ryiov  approxirm  tion  (17.9)  is  limited  by  (17.13),  hence 
it  is  valid  also  in  I  and  Illy.  It  is  somewhat  less  restricted  in 
Illy  than  the  WKB  approximation  to  the  left  part  of  this  region, 

(iii)  The  radio-wave  results  (18.20)  are  restricted  by  L  >>  kL02  and  by 

(18.27)  hence  they  hold  in  II  and  IVg.  Thus,  the  Rice  distribution 
of  I  appears  to  be  restricted  to  radiowave  propagation. 

(iv)  The  coherent  wave  (14.6).  This  result  is  the  most  general  one. 

It  holds  in  all  regions  I -IV,  including  IVgg.  The  point  is  that 
(13.5)  can  be  derived  from  (13.4)  without  using  a  simplified  ex¬ 
ponential  in  (5.4)  because  /k2  =  0  for  odd  m  and  Qm2/k2  =  K^./k2 
for  even  m.  Therefore  no  cumulative  error  here. 

(v)  The  mu tual-oohei-enoe  factor  (1S.S).  This  result,  too,  holds  every¬ 
where  in  I-IV  and  it  is  also  unrestricted  by  the  sagittal  approxi¬ 
mation  because  after  summing  out  beads  (not  restricted  by  sagittal 
errors!)  we  can  repe  the  derivation  of  (15.1)  from  (13.6)  using 
the  full  exponential  of  (5.4). 

Finally,  we  comment  on  the  significance  of  the  two  subregions  of  III. 


The  weak- optics  region  Illy-  This  is  the  region  where  modif ied-Rytov  holds, 


and  where  o( 
variance  vs. 
log-normal. 


is  small,  it  is  well  understood,  and  it  describes  the  trradiauee- 

L-a£2  region.  The  ir radiance  is 


a,/  curves  in  the  rising  small- 


The  strong-oi  tics  region  Ills •  This  is  the  notorious  region  whe,  the  so-called 
saturation  of  irradiance  variance  as  a  function  of  o£2  is  observed  [11,15].  It 
is  characterized  by  aL  >>  1,  and  1  £  (K^L/k)1'^  a62.  In  the  left  pert  of 
region  Illg,  the  irradiance  variance  has  saturated,  and  I  is  observed  to  be  log¬ 
normal.  It  is  not  clear  from  our  formalism  what  takes  place.  To  the  right  - 
as  one  approaches  the  L  =  kL02  line  -  a  transition  to  the  Rice-distribution  be¬ 
havior  of  region  IV  must  take  place,  but  here  also  we  cannot  write  down  quite 
how  this  takes  place.  The  most  promising  approach  in  terms  of  our  formalism 
is  to  take  QnH-i/k  in^(13.1)  and  to  replace  it  by  0  as  given  in  (8.5),  under 
the  assumption  that  0  also  describes  (8.3)  to  good  approximation.  The  develop¬ 
ment  of  this  work  is  currently  under  way.  No  other  approach  has  yet  led  o 
expressions  from  which  the  probability  distribution  of  I  in  this  region  can  be 
predicted.  The  experimental  results  appear  to  predict  the  empirical  relation¬ 
ship, 


[l  -  exp(-4a  2)J  [l  +  8a  £  °’4] 


(19.1) 


where  8  is  an  unknown  function  of  wind  velocity  (?),  microscale,  wavelength, 
and  perhaps  other  parameters.  The  function  8  appears  to  have  a  magnitude 
slightly  larger  than  unity.  However,  the  relationship  is  empirical,  and  only 
one  out  of  a  number  of  possibilities. 
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